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ABSTRACT 

We study the evolution of the ionization state of the intergalactic medium (IGM) 
at the end of the reionization epoch using moderate resolution spectra of a sample of 
nineteen quasars at 5.74 < z em < 6.42 discovered in the Sloan Digital Sky Survey. Three 
methods are used to trace IGM properties: (a) the evolution of the Gunn-Peterson 
(GP) optical depth in the Lya, /?, and 7 transitions; (b) the distribution of lengths 
of dark absorption gaps, and (c) the size of HII regions around luminous quasars. 
Using this large sample, we find that the evolution of the ionization state of the IGM 
accelerated at 2 > 5.7: the GP optical depth evolution changes from Tq® (1 + z ) 4 - 3 
to (1 + *)* n , and the average length of dark gaps with r > 3.5 increases from < 10 
to > 80 comoving Mpc. The dispersion of IGM properties along different lines of sight 
also increases rapidly, implying fluctuations by a factor of ^ 4 in the UV background at 
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z > 6, when the mean free path of UV photons is comparable to the correlation length 
of the star forming galaxies that are thought to have caused reionization. The mean 
length of dark gaps shows the most dramatic increase at z ~ 6, as well as the largest 
line-of-sight variations. We suggest using dark gap statistics as a powerful probe of 
the ionization state of the IGM at yet higher redshift. The sizes of HII regions around 
luminous quasars decrease rapidly towards higher redshift, suggesting that the neutral 
fraction of the IGM has increased by a factor of ^ 10 from z = 5.7 to 6.4, consistent 
with the value derived from the GP optical depth. The mass-averaged neutral fraction 
is 1 — 4% at z ~ 6.2 based on the GP optical depth and HII region size measurements. 
The observations suggest that z ~ 6 is the end of the overlapping stage of reionization, 
and are inconsistent with a mostly neutral IGM at z ~ 6, as indicated by the finite 
length of dark absorption gaps. 


Subject headings: quasars: general; cosmology: observations; quasars: absorption lines; 
intergalactic medium 


1. Introduction 

After the recombination epoch at z ~ 1100, the universe became mostly neutral, until the first 
generation of stars and quasars reionized the intergalactic medium (IGM) and ended the cosmic 
“dark ages” (e.g., Rees 1998). Cosmological models predict reionization at redshifts between 6 
and 20 (e.g., Gnedin & Ostriker 1997, Gnedin 2000, 2004, Ciardi et al. 2001, Benson et al. 2001, 
Razoumov et al. 2002, Cen 2003a,b, Ciardi, Ferrara & White 2003, Hairnan & Holder 2003, 
Wyithe & Loeb 2003a, Somerville et al. 2003, Choudhury & Ferrara 2005). When and how the 
universe reionized remains one of the fundamental questions of modern cosmology (for reviews of 
theoretical models, see Barkana &; Loeb 2001, Loeb &; Barkana 2001, and Ciardi & Ferrara 2005). 

The last few years have witnessed the first direct observational constraints on the history of 
reionization. Lack of complete Gunn-Peterson (GP, 1965) absorption at z < 6 indicates that the 
IGM is highly ionized by that epoch (e.g. Fan et al. 2000, 2001, Becker et al. 2001, Djorgovski et 
al. 2001, Songaila & Cowie 2002). Spectroscopic observations of the highest-redshift quasars and 
galaxies provide a number of observational constraints suggesting that z ~ 6 indeed marks the end 
of the reionization epoch: 

• The detection of complete GP absorption troughs in the spectra of quasars at z > 6 suggests 
that the neutral fraction of the IGM is increasing rapidly with redshift (Cen & McDonald 
2002, Fan et al. 2002, Lidz et al. 2002, White et al. 2003, Gnedin 2004). This feature is 
consistent with the phase of IGM evolution at which individual HII regions overlapped; 
at this point, the IGM ionization state evolved rapidly, similar to a phase transition (cf. 
Songaila 2004 and discussion below). The lower limit of the volume averaged neutral fraction 
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is ~ 10 3 based on Gunn-Peterson measurements, although it could be much higher in the 
deepest troughs. 

• The sizes of large, highly ionized regions around luminous quasars at z > 6 are consistent 
with those predicted for a quasar “Stromgren sphere” expanding into a surrounding largely 
neutral IGM (Mesinger & Haiman 2004, Mesinger et al. 2004, Wyithe & Loeb 2004a). The 
neutral fraction of the IGM is estimated to be as high as > 20% along some lines of sight, 
considerably higher than the lower limit using the Gunn-Peterson optical depth. But these 
measurements are made uncertain because of our lack of knowledge of quasar lifetimes, 
bolometric corrections of quasar luminosities, the bias factor and the dumpiness of the IGM 
around high redshift quasars (e.g., Yu & Lu 2005, Yu 2005). 

• The observed temperature evolution of the IGM at z ~ 3 — 5 suggests that the reionization 
epoch is not much higher than 8 (Hui & Haiman 2003). In models with reionization at much 
higher redshift, the IGM cools to a much lower temperature than that observed at z ~ 2 — 4. 

Quasars and AGNs are unlikely to provide enough UV photons to reionize the universe at z > 6 
(e.g. Fan et al. 2001, Dijkstra et al. 2004, Meiksin 2005). Star forming galaxies at early epochs 
are the most likely candidates to ionize the IGM, although current observations still have poor 
constraints on the total ionizing photon budget at z ~ 6 (Yan & Windhorst 2004b, Stiavelli et al. 

2004) . Recent measurements of the luminosity density of Lyman break galaxies using deep HST 
observations show a moderate decline in total UV luminosity of star forming galaxies at z > 6 
(Bunker et al. 2004, Yan & Windhorst 2004b, Bouwens et al. 2004a,b). If the reionization process 
was similar to a phase transition, the high-redshift quasar observations suggest that the Universe 
could have been mostly neutral as late as z = 6 — 8. 

Polarization measurements of the cosmic microwave background (CMB) based on WMAP 
first-year data (Kogut et al. 2003) show a large optical depth due to Thompson scattering of 
electrons (r = 0.17 ± 0.04) in the early Universe, suggesting that the IGM was largely ionized 
by Z rN -' 17 ± 4. WMAP three-year data (Spergel et al. 2006) show a smaller optical depth 
(r = 0.09 ± 0.03) based on the new EE measurements, indicating a larely ionized IGM by 
z ~ 10 ± 3. The WMAP polarization results, combined with low redshift constraints, suggest that 
reionization could be an extended process rather than a phase transition (Cen 2003a, b, Wyithe & 
Loeb 2003a, Haiman & Holder 2003, Somerville et al. 2003, Gnedin 2004, Choudhury & Ferrara 

2005) . Even at z ~ 6, several lines of evidence suggest that a uniform end of reionization with a 
clean phase transition in the ionization state is too simplistic: 

• While the optical depth no doubt has increased at z > 6, it is not agreed whether the 
increase is similar to a sharp phase transition (e.g., Fan et al. 2002), or is simply a 
continuation of the gradual thickening of the Lycc forest (e.g. Songaila 2004). Although 
optical depth measurements of quasar spectra are consistent among different studies, the 
different interpretations arise because: (1) complete Gunn-Peterson troughs give only a 
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lower limit to the optical depth; ( 2 ) there are complications in the interpretation of the 
Ly/3 measurements; and (3) there is a lack of a clear definition of “gradual” vs. “sharp” 
transitions. 

• There are large line of sight variations from one quasar to another. The GP trough detected 
in the spectrum of SDSS J1030+0524 (z = 6.28, Fan et al. 2001, Becker et al. 2001, White 
et al. 2003) shows complete absorption in the Lya, (3 and 7 transitions. However, the line of 
sight of SDSS J1148+5251 (z = 6.42, Fan et al. 2002) shows clear transmission, especially in 
the Ly (3 and 7 transitions (White et al. 2003, 2005, Oh & Furlanetto 2005), indicating that 
the IGM along this line of sight is still more than 99% ionized. 

• The lack of evolution in the luminosity functions of Lya emitting galaxies from z ~ 5.7 to 
z ~ 6.5 is also consistent with the picture that the IGM is largely ionized at z ~ 6 ; in a 
largely neutral IGM, an extended Gunn-Peterson absorption damping wing would attenuate 
the Lya emission of high-redshift galaxies (Hu et al. 2004, 2005, Malhotra & Rhoads 2004, 
Stern et al. 2005, Malhotra & Rhoads 2006). However, due to uncertainties in the clustering 
of Lya galaxies and the dumpiness of the IGM, a relatively high neutral fraction cannot yet 
be ruled out by the Lya emitter observations (Santos 2004, Haiman & Cen 2005, Wyithe et 
al. 2005, Furlanetto et al. 2006). 

Reionization is a complex process which results from the interplay of structure formation, 
early star formation and feedback, and radiative transfer in a clumpy IGM. Recent developments 
in both semi-analytic models and simulations aim at understanding these processes in the context 
of improved observational constraints and future observations. Models that include feedback from 
the formation of the earliest, metal-free populations can generally fit both the Gunn-Peterson 
measurements at low redshift and a high CMB polarization measurement (e.g., Cen 2003a, b, 
Wyithe & Loeb 2003a, Haiman & Holder 2003, Melchiorri et al. 2005). Cosmological simulations 
with improved resolution, box size and treatment of radiative transfer (e.g. Gnedin 2000, 2004, 
Razoumov et al. 2002, Paschos & Norman 2005) are beginning to make realistic comparison with 
observations possible. 

To make further progress observationally, spectroscopy of a large sample of luminous sources 
is needed to explore the line of sight variation in the GP optical depth measurements expected due 
to dumpiness in the IGM and clustering of ionizing sources (e.g. Wyithe & Loeb 2004b, Paschos 
& Norman 2005). Fan et al. (2002, Paper I) presented measurements of the ionization state of 
the IGM at z ~ 6 and constraints on the epoch of reionization, using a sample of four quasars at 
z = 5.74 — 6.28. In this follow-up paper, we describe measurements of Lyman absorption using 
spectroscopic observations of a sample of nineteen quasars at z = 5.74 — 6.42, discovered from the 
imaging data of the Sloan Digital Sky Survey (Fan et al. 2000, 2001, 2002, 2003, 2004, 2006). We 
address the following three questions: 


1. How fast do the average and the dispersion of the IGM ionization state evolve in the 
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redshift range z = 5 — 6.4? In particular, is there a well-defined break in these quantities as 
a function of redshift? 

2. Is the IGM consistent with a large neutral fraction at least along some lines of sight? What 
is the upper limit of the neutral fraction when considering all the observational constraints? 

3. Can we define a set of statistics from the observed spectrum that are also easily extractable 
from cosmological simulations and allow direct comparison with theoretical models (e.g., 
Songaila Sz Cowie 2002, Paschos & Norman 2005)? 

The paper is organized as follows: in §2, we describe our spectroscopic sample of 19 SDSS 
quasars at z = 5.74 — 6.42. In §3, we measure the evolution of Gunn-Peterson optical depths. We 
present Lya measurements in §3.1, discuss in detail the optimal way of combining Lya, (3 and 
7 measurements in §3.2, and present the measurement of the dispersion of optical depth along 
different lines of sight in §3.3. In §4, we derive a number of quantities characterizing the ionization 
state of the IGM: the evolution of the ionizing background based on a photoionization model 
(§4.1), and of the mean free path of ionizing photons and the neutral fraction of the IGM (§4.2), 
following the methods developed in Paper I. We compare these results with simulations in §4.3. 
In §5, we study the presence of complete Gunn-Peterson absorption troughs in quasar spectra 
(§5.1), and discuss in detail the statistics of the dark gap distribution (§5.2). In §5.3, we use these 
statistics to place an independent upper limit on the IGM neutral fraction. In §6, we measure 
the sizes of HII regions around these 19 quasars, and use their evolution to place an independent 
constraint on the rate of neutral fraction evolution at z ~ 6 . Finally, we summarize the results 
in §7. Through the paper, we use the WMAP cosmology (Ho = 71 km s -1 Mpc^ 1 , Ha = 0.73, 
Om = 0.27, and = 0.04, Spergel et al. 2003) to present our results. 


2. A Spectroscopic Sample of z ~ 6 Quasars 

The Sloan Digital Sky Survey (York et al. 2000) is using a dedicated 2.5m telescope (Gunn et 
al. 2006) and a large format CCD camera (Gunn et al. 1998) at the Apache Point Observatory in 
New Mexico to obtain images in five broad bands (it, g , r, i and z, centered at 3551, 4686, 6166, 
7480 and 8932 A, respectively; Fukugita et al. 1996) of high Galactic latitude sky in the Northern 
Galactic Cap. About 7500 deg 2 of sky in the Northern Galactic Cap have been imaged at the time 
of this writing (Oct 2005). The imaging data are processed by the astrometric pipeline (Pier et 
al. 2003) and the photometric pipeline (Lupton et al. 2001), and are photometrically calibrated to 
a standard star network (Smith et al. 2002; see also Hogg et al. 2001, Ivezic et al. 2004, Tucker et 
al. 2005). 

Quasars with redshift greater than 5.7 are characterized by very red i — z colors. Since 2000, 
we have been using the imaging data from the SDSS to select i-dropout candidates. At z > 5.7, 
all quasars have deep Lyct absorption and satisfy our color-selection criteria (i — z > 2.2); the 
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Fig. 1.— Spectra of our sample of nineteen SDSS quasars at 5.74 < z < 6.42. Twelve of the spectra 
were taken with Keck/ESI, while the others were observed with the MMT/Red Channel and Kitt 
Peak 4-meter/MARS spectrographs. See Table 1 for detailed information. 
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color-selection will not bias IGM measurements. Because such quasars are very red and faint, and 
can be confused with cosmic rays and brown dwarfs, we do not use the automated SDSS quasar 
selection algorithm (Richards et al. 2002a) to select and target them. In a series of papers (Fan 
et al. 2000, 2001, 2003, 2004, 2006), we have reported the selection and identification of nineteen 
quasars selected from about 6600 deg 2 of SDSS imaging. The sample covers the redshift range 
2 i = 5.74 — 6.42, and a range of absolute magnitudes in rest-frame 1450A, M 1450 = —26.2 to —27.9. 
They are the most distant luminous quasars published to date (see also Mahabal et al. 2005, R. 
Cool et al. in preparation). 

Table 1 summarizes the properties and spectroscopic observations of these quasars. These 
spectra were taken with a number of different spectrographs and exposure times, resulting in a 
range of signal-to-noise ratio (S/N) and spectral resolution. Twelve objects in the sample have 
spectra taken with the Echellette Imaging Spectrograph (ESI, Epps & Miller 1998) on the Keck 
II telescope, of which four are presented in this paper for the first time (see White et al. 2003 for 
details of observations and data reduction). These data have spectral resolution R ~ 3000 — 6000, 
depending on the seeing and slit used. To present uniform data, we have binned all the Keck/ESI 
spectra to a resolution of R = 2600. Although the S/N of these spectra vary by a factor of ~ 7, 
we have sufficiently long observations for all quasars at z > 6.1 to ensure proper measurement 
of complete Gunn-Peterson troughs. For quasars at z < 6.1, and absorption redshift z a b s < 5.7, 
the IGM transmission is well detected even with relatively short exposures. In most of this 
paper, we will use averaged transmission and effective optical depth as the main observables. The 
non-uniform S/N among the lower redshift quasars is not a limiting factor in the analysis: in fact, 
uncertainties in the average transmission are dominated by large sample variance (§3.3). Thus 
our data are adequate for our purposes, although calculating higher-order statistics requires more 
uniformly high S/N spectra, such as presented by White et al. (2003), Songaila (2004) and Becker 
et al. (2006). 


3. Evolution of Gunn-Peterson Optical Depth 

Figure 1 presents the spectra of the 19 quasars. At least two of the objects (SDSS J1048+4637 
12 , 2: = 6.20, Maiolino et al. 2004, and SDSS J1044-0125, z = 5.74, Djorgovski et al. 2001, 
Goodrich et al. 2001) are broad absorption line (BAL) quasars. In analyzing their spectra, we have 
excluded regions that are affected by the Lya and Ly (3 BAL features. Our optical spectroscopy 
is usually adequate to detect SilV BAL troughs. However, at 2: ~ 6 , the strongest CIV BAL 
troughs are redshifted to A ^ 1 /jrn. Fewer than half of the objects in the sample have adequate IR 
spectroscopic observations to detect the CIV BAL feature. Note that the fraction of BAL quasars 
at 2: < 4 is ~ 15% in color-selected samples (e.g. Reichard et al. 2003). 


12 Throughout this paper, we use the shortened version of quasar name, SDSS Jhhmin+ddmm. Complete IAU 
names are given in Table 1. 



Figure 1 shows strong redshift evolution of the transmission of the IGM: transmitted flux is 
clearly detected in the spectra of quasars at z < 6 and blueward of the Lya emission line; the 
absorption troughs deepen for the high-redshift quasars, and complete Gunn-Peterson absorption 
begins to appear along lines of sight at z > 6.1. In this section, we present detailed measurements 
of the average Gunn-Peterson optical depth in the Lya, /3 and 7 transitions and discuss the 
variation of the optical depth along different lines of sight. 


3.1. Lya 


The Gunn-Peterson (1965) optical depth to Lya photons is 


tgp =- f a \ a H~ 1 {z)nm, ( 1 ) 

m e c 

where f a is the oscillator strength of the Lya transition, A Q = 1216A, H(z) is the Hubble constant 
at redshift z. and nni is the density of neutral hydrogen in the IGM. At high redshifts, using 
H(z) oc hflm 2 ( 1 + z) 3 / 2 , and converting hydrogen density to the present-day baryon density 
parameter the GP optical depth for a uniformly distributed IGM can be re-written as: 


tgp ( z ) 


1.8 x 10 5 /i _1 H" 1/2 





3/2 



( 2 ) 


We define the transmitted flux ratio as the average ratio of observed flux to that in the absence of 
absorption: 

Tfeb.) = (/;'- U'T), (3) 

where the average is over a certain wavelength or redshift range along the line of sight. We also 
define the effective GP optical depth, 


tgp = -ln(T). 


(4) 


Note that for a clumpy IGM, most of the flux is transmitted through regions with low jzhp while 
high density regions simply saturate, resulting in a much smaller effective optical depth than the 
average of the pixel-to-pixel optical depth, r cff < (r) (for further discussion, see Paper I, Oh & 
Furlanetto 2005, and §4 of this paper). 

We measure the transmitted flux ratio in the Lya transition as a function of redshift for the 
sample of 19 quasars as follows: (1) we define the maximum redshift z “ ax as the maximum Lya 
absorption redshift which is not affected the proximity effect from the quasar itself (§ 6 ). For the 
two known BAL quasars, regions affected by the BAL are excluded based on their SilV and CIV 
absorption redshifts. For non-BAL quasars, we find z em — z“ ax = 0.07 — 0.23, with the difference 
increasing towards lower redshift and higher luminosity (§ 6 ). (2) The minimum absorption redshift 
considered is chosen as 1 + z“ in = (1 + z era ) x 1040A/A q , where 1040 A is the shortest wavelength 
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that is not affected by the Ly/3+OVI emission line. The redshift interval covered by each line of 
sight ranges from 0.76 to 0.98. (3) We choose a redshift interval Az = 0.15 for the calculation of 
the average in Eq. (3). This interval corresponds to ~ 60 Mpc in co-moving distance at z ~ 6, 
larger than any large-scale structure at that redshift (see §4.1). It is also sufficiently wide that the 
measurement error due to photon noise is always small in the absence of complete GP troughs. (4) 
We assume an underlying intrinsic spectrum with a power law continuum f v oc z/ -0 " 5 (Vanden Berk 
et al. 2001) and a double Gaussian Lya +NV emission line. The transmitted flux is measured 
by extrapolating this intrinsic spectrum /), nt to shorter wavelength. When T approaches zero at 
high redshift, the optical depth measurement depends only logarithmically on the exact shape of 
the power law continuum. As discussed in detail in Paper I, the dominant error of the optical 
depth measurement is sample variance in different regions of the IGM. The Lya +NV emission 
line is only important for the region closest to the Lya emission line in the quasar proximity zone 
(§6). For the whole sample, the Lya transmitted flux measurements cover a redshift range of 
4.9 < z a bs < 6.3, with a total length of 14.6 in redshift (~ 6000 comoving Mpc) when adding up 
all lines of sight. 


Table 2 presents the measurements of the Lya transmitted flux ratio. Figure 2 presents the 
effective optical depth for the Lya transition as a function of redshift. When flux is detected at 
less than the 2-a level, a 2-a lower limit on the optical depth is used. At z a b s < 4.9, we plot the 
data from Songaila (2004), where 95% variation due to sample variance is shown for each redshift 
The dashed curve shows the best-fit power-law for z a b s < 5.5: 


= 0.85 x 



(5) 


The results here are consistent with those in previous papers, now using a significantly larger 
sample at z a b s > 5.7. While at z a b s < 5.7, the evolution of optical depth closely follows the 
(1 + z) 4 ' 3 power law, there appears to be an accelerated evolution at z a b s ^ 5.7, where the optical 
depth rapidly exceeds a simple extrapolation from lower redshift. 


In addition to the increase in the average Lya optical depth, there is also an increase in the 
dispersion of this measurement. Dark regions with Tqp > 5.5 occur as early as z a b s ~ 5.7. Deep, 
complete GP absorption troughs are detected in three different lines of sight with various redshifts 
and lengths (§4.2). For example, at z a bs > 5.9, we have eight independent measurements, resulting 
in two upper limits (tqp >6 — 7) and six detections (r^p = 3.9 — 6.5). The smallest optical depth 
measured at z > 6 does not deviate significantly from the low-redshift extrapolation (see e.g. 
Songaila 2004). We discuss the increasing dispersion of optical depth in §3.3. 


3.2. Ly/3 and Lyy 

At a given absorption redshift, the GP optical depth tgp oc /A, where / and A are the 
oscillator strength and rest-frame wavelength of the transition. For the same neutral hydrogen 
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Z abs 

Fig. 2.— Evolution of the Lya Gunn-Peterson optical depth with redshift averaged over intervals 
A z = 0.15 along each line of sight. At 2: a b s = 4.8 — 6.3, the sample of 19 quasars in this paper 
yields a total of 97 independent measurements covering a total redshift interval of = 14.6 
(large symbols). In complete GP troughs in which no flux is detected, the 2 -a lower limit on 
optical depth is indicated with an arrow. We also include the measurements at lower redshift 
from Songaila (2004, small symbols). The dashed curve shows the best-fit power-law for z a i, s < 5.5: 
Tqp = (0.85A0.06) x ( ±±2) . At z a b s ^ 5.7, the evolution accelerates, with increased dispersion 

and a rapid deviation from the power-law relation. 
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density, the GP optical depths of Ly/3 and Lyy are factors of 6.2 and 17.9 smaller than that of 
Lya, respectively. When Lya absorption saturates at z > 6, Ly/3 and Lyy can in principle provide 
more stringent constraints on the IGM ionization state. However, just as we have seen that the 
average and effective GP optical depths are different, in a clumpy IGM, the ratio of effective 
optical depths between Lyct and Ly/3 / Lyy transition are smaller than the factors above. The 
Ly/3 /Lyy absorptions in the quasar spectra are further affected by the lower order transitions 
in the foreground (at lower z a bs) which have to be removed when measuring optical depth. In 
this subsection, we first present direct measurements of Ly/3 absorption, then discuss the ratio of 
optical depths among different transitions. Lyy measurements are discussed at the end of this 
subsection. 

We measure the Ly/3 transmission using the same procedure as for Lya. The minimum 
redshifts for the Ly/3 measurements are chosen as: 1 + = (1 + z em ) x 970A/A p to avoid 

contamination from Lyy absorption. Because of this, 15 of the 19 lines of sight have only one 
independent measurement at A z = 0.15. Two lines of sight have two independent data points, 
while two others have an uncontaminated Ly/3 redshift interval narrower than 0.15. The Ly/3 
measurements cover a redshift range of 5.3 < z a b s < 6.3 with a total redshift interval of 2.9. 
Table 3 presents the observed transmitted flux ratio in the Ly/3 region, without correcting 
for the foreground Lya absorption. Figure 3 presents the GP optical depth due to Ly/3 
absorption along the nineteen lines of sight in our sample, as well as lower redshift measurements 
from Songaila (2004). In Figure 3, we have corrected for the foreground Lya absorption at 
z abs = (1 + Zabs),\g/A a — 1. using Eq. (5) above. Note that while this correction is statistically 
correct, it increases the dispersion of the Ly/3 GP optical depth. The dashed line in Figure 3 is the 
best-fit relation at z a b s < 5.5: 

_i_ y \ 4.3±0.3 

ir) • < 6 » 

We fixed the power-law index to be the same as Eq. (5) in the fit. Compared to Figure 2, the 
deviation from the power-law relation at z a b s > 5.7 is more pronounced in term of the ratio of 
optical depths: at z a b s ~ 6, the lower limit on the average Ly/3 optical depth is ~ 3.5, compared 
to the extrapolated value of 1.7 from Eq. (6). Also note that among the three Lya lower limits 
in Figure 2 that also have valid Ly/3 measurements, one has a clear detection of flux in Ly/3 , 
indicating that Ly/3 is indeed more sensitive to large neutral fraction than Lya. The two other 
Ly/3 points, however, still do not have detectable flux, suggesting a high neutral fraction. 

A number of authors have studied the evolution of the ionizing background and neutral 
fraction of the IGM by combining Lya and Ly/3 measurements (Paper I, Cen Sz McDonald 2002, 
Lidz et al. 2002). They show that the Ly/3 measurements indeed yield stronger constraints on the 
ionizing background at z a b s > 6. Songaila (2004) and Oh & Furlanetto (2005) explicitly estimated 
the ratio of effective optical depths in the two transitions in a clumpy IGM. Here we expand the 
discussion in Paper I, and calculate the ratios using two independent methods. Note that in Paper 
I, the difference between the Lya and Ly/3 wavelengths was neglected in the calculations, resulting 


eff 

T 0 


= (0.38 ± 0.05) x 
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Fig. 3.— Evolution of the Ly/3 Gunn-Peterson optical depth with redshift. At z a b s = 5.3 — 6.3, 
there are 19 independent measurements covering a total redshift interval of A z = 2.9. These 
measurements have been corrected for the foreground Ly a absorption. We also plot the 
measurements at lower redshift from Songaila (2004). The dashed curve shows the best-fit power- 
law for z a -b s < 5.5: rjf = 0.38 x • Like Lya , at z a b s ^ 5.7, the evolution accelerates with 

a rapid deviation from the power-law relation. 
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in an estimated ionizing background that was 7% too high at z a b s > 6 (see §4.1). 


We first estimate the conversion factor using a model in which the IGM is isothermal and 
photo-ionized by an uniform UV background, following Paper I and similar treatments by other 
authors. A detailed discussion of the formalism can be found in Paper I, and will only be 
described briefly here. We define the fractional density of the IGM as A = p/{p). Assuming local 
photoionization equilibrium, whereby the netural hydrogen fraction depends on the square of the 
local density, one can show that for a region of IGM with density A, 


t(A) oc 


(i + z) 4 - 5 qW) 2 a[r(A)] 2 
hT(A,z)n^ 


(7) 


where T is the photoionization rate, assumed to be uniform here, and a(T) is the recombination 
coefficient at temperature T (Abel et al. 1997), a(T ) = 4.2 x 10 _ 13 (T/10 4 K) _ 0 - 7 cm 3 s _1 . The 
relation T(A) oc A 7 is the equation of state of the IGM, where 7 = 0 — 1. We assume 7 = 0 for 
isothermal case here. The observed transmitted flux ratio or effective optical depth is averaged 
over the entire IGM density distribution, 

f OO 

T = {e~ T )= e~ T ^p(A)d(A), ( 8 ) 

Jo 

where p{ A) is the distribution function of the density of the IGM. We calculate T in a clumpy 
IGM using a parametric form for the volume-weighted density distribution function p( A) 
(Miralda-Escude, Haehnelt, & Rees 2000): 


p( A) = A exp 


(A -2 / 3 - C 0 ) 2 

2(2<5 0 /3 ) 2 


A -/3 , 


(9) 


where do = 7.61/(1 + z), and /?, Co and A are numerical constants given in Table 1 of Miralda- 
Escude et al. (2000) at several redshifts, derived from their simulation which assumed a cosmology 
consistent with that of WMAP. We linearly interpolate these constants as a function of redshift. 
Eq. (7) is appropriate for all Lyman series lines, with different values of the proportionality 
constant. Eq. (7) also shows that if the ionizing background and the dumpiness of the IGM do 
not evolve with redshift, then r oc (1 + z) 4 ' 5 . At z a b s < 5.5, this is the leading term in the observed 
scaling relation (Eqs. 5 and 6 ), and the increasing dumpiness and decreasing ionizing background 
as one moves to lower redshift roughly cancel each other out. A significant deviation from this 
relation (§3.3) indicates an accelerated evolution in the ionizing background. 


We calculate the GP optical depth in different Lyman transitions as a function of IGM 
photoionization rate at z = 6. Figure 4 shows the conversion factors and as a 

function of the photoionization rate r_i 2 , in units of 10 “ 12 s _1 . The conversion factor shows a mild 
dependence on the ionizing background. For parameters of interest (see §4.1), r_i 2 = 0.01 — 0.25, 
the Lycc - Ly/3 conversion factor T® ff /rJ ff = 2.5 — 2.9, and the Lya - Lyj conversion factor 
T a / T 7 = 4.4 — 5.7. These values are factors of 2 — 4 smaller than those in the case of a uniform 
IGM, consistent with those found by Oh & Furlanetto (2005). 
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Fig. 4.— The optical depth conversion factors between Lyman series transitions as a function of 
the UV background at z = 6. For a clumpy, isothermal IGM and a uniform background, the Lya 
- Ly/3 conversion factor lies in the range of 2.5 - 2.9, while the Lya - Lyy conversion factor lies 
between 4.4 and 5.7. 
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The calculation above is subject to the assumptions we made about the dumpiness of the 
IGM, its equation of state, and the uniformity of the ionizing background. We also calibrate the 
Lya - Ly (5 conversion factor empirically: we calculate the average transmitted flux ratios (not the 
average optical depth) in Lya and Ly/3 over the same redshift range, T a and Tp, the latter corrected 
for foreground Lya absorption, and define the conversion factor as: t®®/t|® = ln(Tp/T a ). This 
method simply requires that the Lya and Ly/3 measurements yield consistent results on the 
optical depth and neutral fraction. We find that for z a bs = 5.4 — 5.8: t°®/t,!)® = 2.25, and for 
Zabs > 5.9, t®®/tJ® = 2.19. This conversion factor is somewhat smaller than that found using 
the photo-ionization model. Oh & Furlanetto (2005) show that a more clumpy IGM, or a UV 
background in which T is correlated with density, results in a smaller conversion factor. For the 
remainder of this paper, we adopt the conversion factors: = 2.25, and t®®/t®® = 4.4. We 

will use the notation Tq® as the effective optical depth converted to the Lya transition unless 
otherwise noted. 

The measurement of Lyy optical depth is severely restricted by the presence of Ly<5 absorption 
at the low-redshift end and the quasar proximity effect at the high-redshift end. Only the three 
lines of sight towards the three highest redshift quasars at z > 6.2 provide a redshift range 
Ax > 0.06 for this measurement. As before, the Ly 7 optical depth has been corrected for 
foreground absorption in Lya and Ly/3 , using Eqs. (5) and ( 6 ) and converted to the Lya values. 
The results are summarized in Table 4. In all three cases, the Ly/3 regions probed correspond to 
redshift ranges with deep Lya GP troughs. 


3.3. Accelerated Evolution of the Average and Scatter of GP Optical Depth 

Figure 5 shows the combined measurements of optical depth evolution using the Lya, Ly (3 
and Lyy transitions, with Ly/3 and Lyy measurements converted to their Lya equivalent assuming 
conversion factors of 2.25 and 4.4, respectively. For x a b s = 4.9 — 6.3, we have also calculated the 
mean and standard deviations of Tq® among different lines of sight at a redshift interval of 0.2 
based on Lya and Ly/3 measurements. When calculating the average and standard deviation, for 
points where both Lya and Ly/3 measurements are available, we use the Lya measurements in 
the cases where there is detected Lya transmitted flux, and use the Ly/3 measurements or lower 
limits when no Lya flux is detected, regardless whether Ly/3 flux is detected. When there is a 
non-detection, we simply use the 2 -<7 lower limit values to calculate the average and standard 
deviation. Therefore, they are lower limits for the two highest redshift bins. The results are 
tabulated in Table 5. We do not use Lyy measurements in Figure 5 to calculate the average and 
standard deviation, because they only cover a small redshift range. Including Lyy measurements 
increases the values of both the average and the scatter of the GP optical depth. 

From the figure, we note that: 

(1) At x < 6, the Lya and Ly/3 measurements are fully consistent. This is a result of the 
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Fig. 5.— Evolution of optical depth with combined the Ly a, Ly/3 and Ly 7 results. The Ly/3 
measurements are converted to Lya GP optical depth using a conversion factor of 2.25. The values 
in the two highest redshift bins are lower limits since they both contain complete GP troughs. 
The dashed line is a redshift evolution of oc (1 + z) 4,3 . At z > 5.5, the best fit evolution has 
Tqp oc (1 + c:) >10 ' 9 , indicating an accelerated evolution. The large filled symbols with error bars 
are the average and standard deviation of optical depth at each redshift. The sample variance also 
increases rapidly with redshift. We also include the Lyy measurements (square symbol, Table 4); 
they cover only the high redshift end of the Lya and Ly/3 bins due to contamination from Ly<5 
absorption along the line of sight at lower redshift. 
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calibration procedure used above. At the highest redshift when no Lyct flux is detected, the Ly/3 
optical depth measurement either gives a secure optical depth measurement that is somewhat 
higher than the Lya lower limit (in one case), or provides a much more stringent lower limit (in 
two cases). Note that the Lyy measurements only cover the high redshift end of each bin at which 
the absorption is particularly strong. 

(2) With the additional constraints from Ly/3 measurements, the trend of accelerated evolution 
in optical depth is clearly established: at z a b s = 5.85, the lower limit on the averaged optical 
depths is 4.0, compared to the extrapolated value of 3.3 based on Eq. (5); and at z a b s = 6.1, the 
lower limit on the averaged optical depths is 7.1, compared to the extrapolated value of 3.8. We 
fit a power-law for redshifts z a bs >5.5, and find that 

r GP a (l + z) >10 ' 9 , Aibs = 5.5 - 6.3. (10) 

This should be compared with tgp oc (1 + z) 4,3 at lower redshift, indicating a rapid change in the 
evolution of optical depth. 

(3) The accelerated evolution of the optical depth is accompanied by an increase in the 
dispersion of the optical depth: at z a b s < 5.0, the standard deviation a(r ) ~ 0.3. It increases to 
c t(t) ~ 0.6 at Zabs = 5.5, and ct(t) ~ 0.8 at z a bs = 5.8. At z a bs > 6.0, the lower limit on u(r) is 
2.1. The increased dispersion z a b s > 5.7 is due to the fact that some lines of sight have complete 
GP troughs, with others have optical depths only mildly higher than the power-law extrapolation 
from low redshift. Different data points are independent measurements of the optical depth in 
different IGM regions (§4.1). This result suggests that the ionization state of the IGM has large 
sample variance at the end of reionization (§5). 

Using a log-normal density distribution, and carrying out the integral of Eqs. (7), (8), and 
(9), Songaila & Cowie (2002) found: 

T& = AW+WV - h (A) + const, (11) 

where A oc 1/T, and (3 = 2 — OTy for an uniform radiation background (see also Oh Sz Furlanetto 
2005). If the IGM is isothermal (/3 = 2), we have r cff oc T~ 0 ' 4 as the leading term. This scaling 
is not sensitive to the detailed IGM density distribution (Eq. 9). Therefore, the evolution of the 
GP optical depth is driven by the rate of evolution in the UV background: if we observe that the 
GP optical depth increases by ~ 2.5 from z ~ 5.5 to 6.2, we can infer that the UV background 
decreases roughly by an order of magnitude. We carry out a detailed calculation in the next 
section. 


4. Evolution of the Ionization State of the IGM 

In Paper I, we used the evolution of the GP optical depth and the formalism described 
in Eqs. (7) - (9) to calculate a number of IGM ionization parameters, including the ionizing 
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background and neutral fraction of the IGM. We also calculated the mean free path of ionizing 
photons using the model of Miralda-Escude, Haehnelt, & Rees (2000). Here we carry out similar 
calculations with the new dataset, and introduce an alternative way of deriving the mean free 
path. However, we caution that the ionization parameters derived in this section are model 
dependent: they assume a particular model of the IGM density distribution and photoionization 
equilibrium with a uniform ionizing background. We will discuss the limitations of these derived 
properties in detail in subsequent subsections. 


4.1. Ionizing Background 

Figure 6 shows the estimated photoionization rate evolution with redshift based on Eqs. (7) 
to (9) and the GP optical depth measurements in Figure 5. Details of these calculations are 
discussed by McDonald & Miralda-Escude (2001) and in Paper I. Consistent with the trend seen 
in Paper I, the drop of ionizing background accelerates at z a b s > 5.7: from z a b s = 5.6 to 6.1, V 
declines by a factor of > 4. In the deepest GP troughs, the 2-cr upper limit of the photoionizing 
rate derived from the Ly (3 transition is T-i 2 < 0.025. The most stringent constraint comes from 
the Lyy measurement (Table 4) in quasar J1030+0524 at z a bs = 6.11 — 6.17, where tqp > 14.2 
(2-i t) is consistent with a photoionizing rate r_i 2 < 0.012 in that region of the IGM, more than 
an order of magnitude lower than the average background at z a b s = 5.5. Also note the relatively 
flat evolution between z a b s = 5 and 5.5. This lack of evolution in T was first noted by Cen & 
McDonald (2002) in a small sample of SDSS quasars. The IGM ionizing background is expected 
to increase towards lower redshift as the emissivity from star forming galaxies and AGNs increases 
and the photon mean free path increases. Cen & McDonald (2002) suggested that the flattening 
might be due to the change in the star formation rate as the temperature and Jeans mass rises 
shortly after reionization. 

Just as we saw in Figure 5, there is also a rapid increase in the dispersion of the estimates 
of the ionizing background. At x ~ 5, the relative dispersion <r(r)/r is ~ 30%. It increases to 
~ 50% at x ~ 5.5. At x > 6, cr(r)/T > 100% with r_i 2 measurements ranging from < 0.02 to 
~ 0.09, reflecting the fact that some lines of sight have complete GP troughs while others still 
have detectable transmission. In the photoionization model, the observed fluctuation could be due 
to two effects: (1) large scale fluctuations in the underlying density field due to the finite length 
of the line of sight used in the measurement, or (2) intrinsic fluctuations in the UV background 
and/or temperature of the IGM. 

Data at x a b s > 4.8 in Figure 6 are derived from an effective optical depth over Ax = 0.15 
(§4.1), corresponding to 4411” 1 comoving Mpc at x ~ 6. The one-dimensional rms fluctuation of 
the average mass density along a line of sight with length L is given by: 

1 r°° 

<{L) = ^J o dk(F*)W 2 (kR), (12) 
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Fig. 6.— Evolution of the photoionization rate (in units of 10 -12 s -1 ) with redshift. The points 
at z a bs < 4.8 are taken from McDonald & Miralda-Escude (2001). For the highest redshift points, 
measurements using Lya (closed symbols) and Ly j3 (open symbols) transitions are shown separately. 
Large, solid symbols with error bars are the average and 1 -a dispersion of T at high redshift. As 
in Figure 5, we show the Lyy measurements, but do not include them in the calculation of average 
and dispersion in the highest redshift bin. 
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where W(x) = sin x jx is the top-hat window function, and (F?) is the one-dimensional power 
spectrum, related to the normal three-dimensional linear power spectrum as (Kaiser & Peacock 
1991): 

1 poo 

( F i) = -^j k P(y)ydy. (13) 

Using the linear power spectrum determined from the Ly a forest at z ~ 3 by McDonald et 
al. (2005), and assuming that it grows linearly with scale factor at high redshift, we find that 
at z ~ 6, the one-dimensional r.m.s. density fluctuation is a m (L) = 0.11 over A z = 0.15. From 
Eq. (7), r oc pi for a given optical depth; a 10% fluctuation in the average density leads to a 
fluctuation of ~ 20% in the resulting T value. This value is somewhat smaller than the dispersion 
of the r measurements at z ~ 5, and is much smaller than that at higher redshift. Over such a 
narrow redshift range, the average IGM temperature is not likely to have changed by more than 
a factor of two (Hui & Hairnan 2003), resulting in a difference in estimated T by < 60%, since 
r oc a(T)/T oc T _a7 /r (§3.2). Therefore, we interpret the large dispersion of optical depth and 
derived ionizing background at z > 5.7 as a result of fluctuations in the UV background near the 
end of reionization: even over the scale size of 44 h _1 comoving Mpc, the average UV background 
fluctuates by a factor of ^ 4. The dispersion of estimated T values at z < 4 is much smaller and is 
mostly due to sample variance of the average density along the line of sight. 

Large fluctuations in the UV background also imply that the transmission is correlated over 
large scales. We calculate the line of sight dispersion of transmitted flux by comparing the average 
transmitted flux over neighboring bins along the same line of sight, and subtract from it the 
difference due to redshift evolution (Eq. 5). We find that at z ^ 5.5, when averaging over redshift 
interval As: < 0.1, the dispersion of line of sight of GP optical depth is ~ 20% smaller than the 
dispersion when comparing different lines of sight. Thus the GP optical depth is correlated over 
scales of tens of comoving Mpc. At A;? > 0.15, this difference is < 10%. At larger A z, the test 
becomes difficult due to the limited line of sight towards any particular quasar. Because of the 
limited sample size, non-uniform S/N of the sample spectra, and strong redshift evolution, we 
do not attempt to calculate the transmission two-point correlation function in this paper (e.g. 
McDonald et al. 2000). By observing high-redshift quasar pairs that have projected separations of 
tens of comoving Mpc, one could check this correlation in the transverse direction (e.g. Mahabal 
et al. 2005). 

Large fluctuations in the UV background are a generic feature of the IGM close to the end of 
reionization, when the individual HII regions start to percolate and most of the UV photons come 
from the earliest highly clustered star forming galaxies. Note that in our photoionization model, 
we assumed a uniform UV background. Here we show that at the end of reionization, there is a 
large background fluctuation. The derived UV background fluctuations will provide a sensitive 
test of models of the reionization history when compared directly with simulations. 
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4.2. Neutral Fraction and Mean Free Path of the IGM 

We use the GP optical depth measured in §3 and the IGM density distribution, Eq. (9), to 
calculate both the volume- and mass-averaged neutral hydrogen fraction. Figure 7 shows the 
estimated volume-averaged neutral fraction using our sample of nineteen quasars. It is compared 
with a recent simulation by Gnedin (2004), which has a sharp transition in the ionization state at 
z ~ 5.5 — 6.5, and a small total Thompson scattering optical depth to CMB r e = 0.06 (compare to 
the WMAP three-year value of r e = 0.09 ± 0.03, Spergel et al. 2006, see §4.3 for more discussion 
on the model comparisons). The strong evolution of neutral fraction in Paper I is confirmed here: 
there is an order of magnitude increase between z ~ 5.1 and 6.1, with the evolution accelerated at 
z > 5.7. We emphasize that the calculation of volume- or mass-averaged neutral fractions depends 
on the density distribution function of the IGM as well as on the assumption of a uniform UV 
background. Oh & Furlanetto (2005) show that most of the IGM transmission is through the most 
highly ionized, underdense regions (see also Figure 7 of Paper I). As discussed in Paper I, the 
volume averaged values are less affected by the ionization state of high-density regions than are 
the mass averaged values. Over a narrow redshift range where the density distribution of the IGM 
evolves little, the ratio of neutral fractions at different redshifts is a robust quantity, especially 
when compared with simulations. 

In Paper I, we calculated the evolution of the mean free path of ionizing photons, using the 
model of Miralda-Escude, Haehnelt, & Rees (2000). They model reionization using the evolution 
of critical overdensity A,;, above which the IGM is still mostly neutral, and below which it is 
mostly ionized. A,; increases as the reionization front progresses to higher density regions. They 
also suggest an expression for the ionizing photon mean free path: A* = Ao[l — -F v (Aj)]~ 2 / 3 , where 
F v (Aj) is the fraction of volume with A < A*, and Ao H = 60 km s^ 1 , which assumes that the 
shapes of isolated density contours remain roughly constant as A,; decreases.. However, the critical 
overdensity A* is not well defined in the simple photoionization picture; we assumed that 95% 
of the total neutral gas lies above Aj in Paper I. In addition, overdense regions in the IGM are 
likely to contain the earliest star-forming galaxies that provide most of the ionizing photons, thus 
they might have been ionized earlier. For these reasons, in this paper, we propose an alternative 
approach to calculate the mean free path A by simply using the mass averaged neutral hydrogen 
density: 

1 / v /o A '' »m(A)p(A)r/A 

A ( v) / 0 %(A)dA ’ ( } 

where (a u ) is the frequency averaged cross section of Lyman limit absorption, assuming an ionizing 
spectrum J v oc z/ -5 (Barkana Sz Loeb 2001), and nni(A) is the neutral fraction corresponding to a 
certain overdensity calculated using Eqs. (7) - (9). The integrals over the density distribution are 
only carried out to A c = 150, above which the region is already virialized and is no longer part of 
the IGM. This expression is accurate in cases in which the mean free path is larger than the scale 
of IGM clustering, or the UV photon field follows the IGM density distribution; these assumptions 
might not be valid at the early stage of reionization. Furthermore, these calculations depend 
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Fig. 7.— Evolution of volume-averaged neutral hydrogen fraction of the IGM. The small dots are 
measurements based on the nineteen high-redshift quasars, while the large points with error bars are 
the means in bins of redshift. The dashed and dotted-dashed lines are the volume-averaged results 
from the simulation of Gnedin (2004) with box sizes of 4 h" -1 Mpc and 8 h^ 1 Mpc, respectively. 
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Fig. 8.— Evolution of the mean free path of ionizing photons inferred from the observations using 
Eq. (13) (solid symbols with error bars) and using the method presented in Paper 1 (dotted line), 
compared with simulations of Gnedin (2004, dashed and dashed-dotted-dashed lines). The mean 
free path at z ~ 5.7 is comparable to or smaller than the correlation length of high-redshift star 
forming galaxies. 
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on the derived neutral fraction and assumed IGM density distribution. While the exact values 
of the mean free path are model-dependent, the strong evolutionary trend is clearly detected. 
Figure 8 shows the mean free path calculated using Eq. (14), where both the mean and 1-a error 
bars are presented. The dotted line shows the mean free path calculated using the method in 
Miralda-Escude et al. (2000) and in Paper I. The two results are fully consistent. They again show 
the accelerated transition at z a b s > 5.7, where only an upper limit can be derived for complete GP 
trough regions. 

At z > 5.7, the mean free path is smaller than 0.5 h _1 proper Mpc (3.5 h _1 comoving Mpc). 
In comoving units, this is comparable to the correlation length of galaxy clustering at z ~ 0, 

(ro ~ 5 h _1 Mpc, Zehavi et al. 2005), and of Lyman break galaxies at z ~ 2 — 4 (ro ~ 4 h _1 
Mpc, Adelberger et al. 2005), and at z = 4 — 5 (ro ~ 4 — 8 h _1 Mpc, Kashikawa et al. 2006). At 
z > 5, evidence of galaxy clustering on scales of ~ 10 Mpc are reported in Hu et al. (2005), and 
Malhotra et al. (2005). High-redshift star forming galaxies, which likely provide most of the UV 
photons for reionization, are highly biased and clustered at physical scales similar to the mean free 
path. The assumption of uniform UV background under which we have derived these results is no 
longer valid. We saw clear evidence for non-uniform T in §4.1. 


4.3. Comparison with Hydrodynamical Simulations 

Hydrodynamical simulations (Gnedin 2002, 2004, Razoumov et al. 2002, Paschos & Norman 
2005) have been used to calculate the evolution of the GP optical depth during the epoch of 
reionization. Gnedin (2004) presented a SLH simulation of reionization by stellar sources with a 
new recipe for radiative transfer using an optically thin variable Eddington tensor approximation 
and box sizes of 4 - 8 h -1 comoving Mpc. In these simulations, the mean mass-averaged neutral 
fraction reaches 10~ 2 and 0.5 at z = 6.3 and 7.3, respectively. The reionization process is modelled 
as similar to a phase transition, completed within A z ~ 1. Similarly, Paschos & Norman (2005) 
presented a set of reionization simulations using the Eulerian code, Enzo, and a box size of 6.8 
h _1 co-moving Mpc. Their default model has a sharp reionization at z re ion ~ 7, around which the 
mass-average neutral fraction increases from 10^ 3 to of order unity. Both simulations have been 
calibrated according to the average GP optical depth measurements at z < 6 using earlier datasets 
(e.g. Becker et al. 2001, Songaila 2004), so the consistency in measurements in this redshift range 
is by design. 

However, these simulations which model the reionization process as a phase transition could 
be too simplistic for two reasons. First, both simulations include only a simple prescription of 
the evolution of emissivity of ionizing sources. The simulations have not captured all the relevant 
physics, especially the effects of early galactic feedback and large scale structure. The total 
Thompson scattering optical depth (r e ) is less than 0.10 in both simulations, inconsistent with the 
WMAP first-year measurements at face value (r = 0.17 ± 0.06, Kogut et al. 2003, Tegmark et al. 
2004), but is consistent with the WMAP three-year measurements (Spergel et al. 2006). A number 
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Fig. 9.— Comparison of the evolution of observed GP optical depth with simulations. Data points 
are the same as in Figure 5. The dotted line is the power-law fit to optical depth evolution at 
z < 5.5: tqp oc (1 + z) 43 . The solid line with error bars is the result from simulations of Paschos 
& Norman (2005), which have 2: re ion ~ 7. The dashed and dotted-dashed lines are the 4 Mpc and 
8 Mpc simulations of Gnedin (2004), which have 2 re ion ~ 6.5. The data points fall between the two 
simulations, but are somewhat closer to those of Gnedin. 
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of semi-analytical analyses (e.g., Cen 2003a,b, Chiu et al. 2003, Wyithe & Loeb 2003a, Haiman 
& Holder 2003 etc.), show that adding feedback from an early generation of galaxy formation, or 
adjusting the star formation efficiency at z > 6 — 10, will produce a prolonged, or multi-epoch 
reionization history. These models predict a GP optical depth consistent with quasar observations 
at z < 6.5, while having a large total Thompson scattering optical depth. 

The second problem is the impact of large scale structure. Both simulations have a limited 
box size (< 10 comoving Mpc). However, Wyithe & Loeb (2004) and Fuiianetto & Oh (2005) have 
studied growth of ionized bubbles at the end of reionization, and estimated that the bubble have 
sizes of > 10 Mpc at the end of the overlapping phase. Wyithe & Loeb (2004, 2005) showed that 
the observed redshift of overlap along different lines-of-sight has a variation of Az ~ 0.15, and the 
GP optical depth variations should be of order unity on a scale of ~ 100 comoving Mpc. These 
results are consistent with the large scale variation of the UV background discussed in §4.1. Thus 
to fully understand the reionization process, simulations with box sizes of the order 100 Mpc are 
needed (e.g., Kohler et al. 2005). 

Figure 9 compares the GP optical depth measurements of Figure 5 with the simulations of 
Gnedin (2004) and of Paschos & Norman (2005). We also compare the estimated IGM neutral 
fraction and mean free path of UV photons with those of the Gnedin simulations in Figures 7 
and 8. In the simulations, the onset of reionization is characterized by an acceleration in the 
evolution of the GP optical depth. Paschos & Norman (2005) showed that when z re ; on — z > 1, the 
GP optical depth scales as a simple power-law tqp oc (1 + z) 4 ' 2 , consistent with what we find at 
z a bs < 5.7. For the last Az = 0.5 — 1, as the epoch of reionization approaches, the optical depth 
evolution accelerates, again mirroring the observed trend. The observed points fall in the range 
covered by the three simulations. A similar rapid change is noted in the volume-averaged neutral 
fraction (Figure 7) and photon mean free path (Figure 8), in which the simulations closely match 
the values estimated from observations. 

The simulations also show that the end of reionization is accompanied by a rapid increase 
of dispersion in optical depth. Paschos & Norman (2005) calculated the flux variance among 
different lines of sight by observing through random lines of sight in the simulation box. The solid 
error bars in Figure 9 are the 68% range of the optical depth measurements at each redshift bin. 
At z < 5.5, the dispersion in the simulation matches the observations well. It also reproduces 
the trend of increased scatter at z > 6, although at somewhat higher redshift. This accelerated 
evolution corresponds to the stage of reionization where individual HII regions begin to overlap 
and the IGM ionization experienced a rapid transition (e.g. Gnedin 2002). The comparison 
of the observations with simulations strongly suggest that z ~ 6 corresponds to the end of the 
overlapping stage of reionization, but for reasons mentioned earlier, it does not provide detailed 
guidance to what happened at much higher redshift, when the reionization history may have been 
more complex. 
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5. Distribution of Dark Gaps 

In the last three sections, we have used the evolution of the effective optical depth to probe 
the ionization state of the IGM. The effective optical depth directly probes the neutral fraction 
and UV ionizing background in photoionization models. However, it is a low order statistic. At 
z > 5, the appearance of the Lya forest is the reverse of that at low redshift: most regions along 
the line of sight are optically thick, while a small fraction of transparent “spikes” provide most 
of the Lya transmission. Meanwhile, ever longer dark absorption gaps - continuous regions with 
large optical depth - emerge at z a h s >5.7. A successful theory of IGM evolution at the end of 
reionization needs to reproduce not only the evolution of the mean optical depth, but also the 
geometry and topology of IGM transmission and the dark gaps. In Pentericci et al. (2002) and 
in Paper I, we examined the distribution of transparent spikes along a few lines of sight. In this 
section, we expand these previous studies and discuss three issues: the presence of complete GP 
absorption troughs at the highest redshift (§5.1); the evolution of the dark gap distribution (§5.2); 
and the use of dark gap lengths to constrain the IGM neutral fraction (§5.3). 


5.1. Gunn-Peterson Troughs in the Highest-Redshift Quasars 

The presence of complete Gunn-Peterson troughs is the most compelling evidence that the 
IGM is approaching the end of the reionization epoch at z ~ 6. To probe the evolution of the IGM 
at z > 6 requires quasars at z > 6.1 to not be affected by the quasar proximity effect. Figure 11 
shows the regions immediately blueward of the Lya and Ly (5 emission lines for the five quasars in 
our sample at z > 6.1; deep absorption troughs with Tq® > 5 and length Az > 0.1 are indicated. 
The ends of troughs are chosen as locations at which significant transmission spikes with r < 2.5 
are detected. Table 6 summarizes the GP optical depth measurements for these troughs. We now 
briefly discuss the case for each quasar: 

SDSS J1148+5251 (z = 6.42). This object is discussed in detail in White et al. (2003, 2005) and 
Oh & Furlanetto (2005). Extended dark regions in Lya absorption start at z a b s > 5.91; however, 
the spectrum shows clear transmissions in the Lya - Lyy transitions at this redshift. White et 
al. (2003) speculated that the Ly (3 transmission could be due to Lya emission from intervening 
galaxies at z ~ 4.9, as suggested by strong CIV absorption at that redshift. Subsequent HST 
imaging indicated that the flux in the Ly/3 trough is a point source coincident with the quasar 
position, inconsistent with the intervening galaxy hypothesis. The Lya and Ly (3 troughs yield 
consistent optical depth measurements: Tq® ~ 5.5 — 7. The consistency between the Lya and 
Ly/3 GP optical depth measurements is direct evidence that at z ~ 6, the IGM is already highly 
ionized at least along this line of sight, although the neutral fraction has increased by a factor of 
~ 4 from that at z ~ 5.7 even in this most transparent line of sight among z > 6.1 quasars. Table 
4 also shows the Lyy measurements in the highest-redshift interval at z = 6.25 — 6.32, where the 
transmission is detected at 3-a level, yielding an effective optical depth: Tq® ~ 15. No Lya or Ly (3 
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Fig. 10. Deep GP troughs in the five highest-redshift quasars known at z > 6.1, in both the Lya 
and Ly/3 transitions. The grey band around zero transmission represents the 1 -a uncertainty of 
the data. The ends of troughs are chosen as locations at which significant transmission spikes are 
detected. SDSS J1048+4637 is a BAL quasar, in which we use the redshift of the CIY BAL feature 
as the upper bound of the GP trough. These quasars show significant line of sight variation in the 
GP optical depth measurements and the redshift where deep troughs first appear. 
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transmission is detectable over this narrow redshift range, consistent with the Ly 7 measurement, 
and suggesting large fluctuations in the neutral fraction along the line of sight. Note that there is 
also a long Lya trough at 5.91 < z a b s < 6.08, although it has clear transmission in Ly (3 , indicating 
a high level of ionization. 

SDSS J1030+0524 ( z = 6.28). This is the first quasar at z > 6 in which a complete GP trough 
was detected (Becker et al. 2001); it still has the longest and deepest absorption trough among 
known quasars, in which no flux in any of the Lyman transitions is detected, yielding the most 
stringent constraint on the evolution of the neutral fraction. 

SDSS J1623+3112 (z = 6.22). This line of sight is very similar to that of SDSS J1030+0524, 
with complete GP troughs in the Lya and Ly (3 transition immediately blueward of the emission 
lines. The starting redshift and length of the complete GP trough are also similar. This object is 
fainter and has a shorter total exposure time, yielding an upper limit 15% weaker than that of 
SDSS J1030+0524. It also has a Lyq trough with no detectable transmission. 

SDSS J1048+4637 (z = 6.20). This is a BAL quasar (Fan et al. 2003, Maiolino et al. 2004). 
The SilV absorption extends to z a bs ~ 5.83, and a somewhat weaker CIV trough reaches to 
^abs = 5.75. We therefore can only measure the IGM properties at lower redshift. The quasar 
has a dark trough at z = 5.63 — 5.75, where no flux is detected in either Lya or Ly (3 at a level 
comparable to SDSS J1623+3112. This is at a redshift considerably lower than the other lines of 
sight. We regard it as less reliable, as the trough could be affected by weaker BAL features as well 
as by absorption troughs from other ions (e.g. PVA1121, Arav et al. 2001). 

SDSS J1250+3130 (z = 6.13). This object shows a dark Lya GP trough at z = 5.69 — 5.85 
with an optical depth Tq® > 4.6 (2-er). The trough doesn’t start immediately blueward of the Lya 
emission line, and has detectable Ly/3 transmission, with Tq® ~ 6.5. The Ly/? measurement is 
comparable to that of SDSS J1148+5251 at considerably higher redshift, although the S/N of the 
spectrum is lower for this quasar. 

The properties of these absorption troughs show significant variation between lines of sight. 
Two of the three quasars at z > 6.2 show complete Lya - Lyy troughs, while the third one has 
only a somewhat higher optical depth than seen at lower redshift. On the other hand, some dark 
troughs already start to appear at z ~ 5.7. We quantify the evolution of the dark absorption gaps 
in the next subsection. 


5.2. Dark Gap Distribution 

Statistics of dark gaps (e.g. Croft 1998) were first applied to a sample of z = 4 — 6 quasars 
by Songaila & Cowie (2002). They studied the evolution of continuous absorption gaps with large 
optical depth (r > 2.5) and found that extended absorption troughs such as those observed by 
Becker al. (2001) in SDSS J1030+0524 can be explained by large variations in the UV background 
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at z > 6. Barkana (2002) presented a semi-analytic model for the length of dark gaps as a function 
of IGM HII region filling factor. Furlanetto, Hernquist, & Zaldarriaga (2004) also presented 
a model of dark gap distribution, considering source clustering and the inhomogeneous IGM. 
Paschos & Norman (2005) and Kohler et al. (2005) studied the dark gap distribution using spectra 
derived from their hydrodynamical simulations. They suggest using the distributions of dark gaps 
as a statistic to compare with observations, and found that the onset of reionization is associated 
with a dramatic increase both in the average and dispersion of gap length. 

Figure 11 shows the distribution of Lya dark gap lengths in the 12 lines of sight that have 
Keck/ESI spectra. We define dark gaps as continuous regions in the quasar spectra in which all 
pixels have optical depth larger than a threshold r m i n . In the upper panel of Figure 11, we use 
T m i n = 2.5, the same as that of Songaila & Cowie (2002) and Paschos & Norman (2005). Both 
the average gap length and its dispersion increase dramatically at z > 5.7. At this optical depth 
threshold, dark regions with Ar ;> 40 comoving Mpc already begin to appear at z a bs ~ 5.3 (e.g. 
Djorgovski et al. 2001). Long gaps with Ar ;> 80 comoving Mpc appear at z a b s ~ 5.7; however, 
they are relatively rare while most lines of sight still have plenty of pixels with r < 2.5 and only 
short dark gaps. At z a b s > 5.9, all lines of sight have long dark gaps with r > 2.5. In many 
cases, the gap length is only a lower limit, as the gap extends all the way to the quasar proximity 
zone. Note that at the high redshift end, imperfect sky subtraction in regions with strong OH 
lines results in large residuals in the spectrum and an underestimate of the dark gap length. This 
transition from isolated, short gaps occupying only a small fraction of pixels, to deep, long troughs 
is the most abrupt transition among all the IGM properties we have presented in this paper, 
and is difficult to explain by a gradual thickening of the Lya forest. Rather, it suggests a phase 
transition in the IGM geometry. Detailed comparison with state-of-art simulations will constrain 
the IGM geometry at this critical stage of reionization. 

Figure 11 compares the observations with the simulations of Paschos Sz Norman (2005) 
described in §4.3. The dashed and dotted lines are model predictions of the average gap length for 
spectral resolution of R = 5300 and R = 36000. Our spectra are binned to R = 2600, and should 
be compared to their low-resolution results. The observations and simulations agree very well at 
z < 5.5. Both also show a strong break in the evolution of dark gap length; however, the break in 
the observations is both earlier (by A z ~ 0.5) and sharper. This is consistent with the comparison 
in Figure 9, and suggests that the reionization epoch is at somewhat lower redshift than in the 
simulations. Note, however, the simulation of Paschos Sz Norman (2005) has a box size of only 6.8 
h -1 comoving Mpc. They used multiple passes through the box to measure the gap sizes larger 
than the box size, which could bias the result. 

At Tmin = 2.5 , all gaps at z > 6.0 are longer than 100 comoving Mpc, and most are only lower 
limits. In order to further constrain the neutral fraction, we use a higher optical depth threshold 
r > 3.5 (shown in the lower panel of Figure 11). Even at z > 6.0, the average length of dark gaps 
is still finite: A?’ ;> 80 comoving Mpc. At this higher threshold, it is clear that the IGM is not yet 
neutral at z 0 # 
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Fig. 11.— Distributions of dark gap lengths, defined as continuous regions in the spectra in which 
all pixels have observed optical depth larger than 2.5 (upper panel) and 3.5 (lower panel) for the 
Ly a transition in the twelve lines of sight that have Keck/ESI observations. Upward arrows are 
gaps immediately blueward of the quasar proximity zone, therefore the length is only a lower limit. 
Solid lines with error bars are mean depth lengths, with 1 -a dispersion given at each redshift bin. 
Long dark gaps start to appear at z ^ 5.6, and the average gap length increases rapidly at z > 5.9. 
In the upper panel, we also compare the observations with simulations of Paschos &; Norman (2005), 
which has a sharp reionization at z ~ 7. The dashed and dotted lines are for spectral resolution 
of 5300 and 30000, respectively, in the simulated spectra. The observed spectra are binned to 
R ~ 2600. The observations appear to have reached the overlap stage at z ~ 6. 
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In the redshift range considered in this paper, the measurement of the dark gap distribution 
and the GP optical depth place similar constraints on the ionization of the IGM. However, at 
z a bs > 6.5, it will become increasingly difficult to map the reionization history using the GP 
optical depth. Currently, with ~ 10 hours of observing time on 10-meter class telescopes, we are 
able to place lower limits on the GP optical depths of r ^ 7 using Lya, and a factor of ~ 2 better 
using Ly (5 and Lyy, corresponding to a neutral fraction of 10~ 2 — 10” 3 . It is difficult to improve 
this measurement further in practice, because the optical depth limit depends logarithmically on 
the observing time and the ability to control sky subtraction systematics. At z > 6, most of the 
pixels are completely dark, and contain little information. After entering the redshift range where 
continuous dark gaps dominate, all useful information is contained in the flux and distribution 
of few remaining transmitting pixels. Therefore, the dark gap distribution provides a powerful 
alternative to map the IGM ionization at high neutral fraction. Paschos & Norman (2005) show 
that even when the average neutral fraction is of the order few tenths, the average gap length is 
still finite, due to the presence of regions that had already been ionized by star forming galaxies. 
In the next subsection, we place an upper limit on the neutral fraction using the finite length of 
the gaps. 


5.3. Constraining the Neutral Fraction with Gap Statistics 


So far in this paper, we have only considered a photoionized IGM without explicitly taking 
into account the presence of strong UV ionizing sources that caused reionization: early star 
forming galaxies and AGNs. These objects create expanding HII regions even in a largely neutral 
IGM, and result in transmitting spikes in the absorption spectrum. However, in a neutral IGM, 
the optical depth at the Lya resonance line center is so large that it generates an extended GP 
damping wing (Miralda-Escude 1998). Small HII regions are still optically thick to the observer 
due to absorption from this damping wing; and only large HII regions around long-lived, luminous 
sources can survive the damping wing absorption. Therefore, the existence of transmitting regions 
between dark gaps, together with a model for the ionizing sources, can be used to place an upper 
limit on the IGM neutral fraction independent of other IGM optical depth constraints. This 
is similar to the constraint using the evolution of Lycc emitting galaxies (Malhotra 2005, Stern 
et al. 2005, Hainran & Cen 2005, Santos 2005). In the case of Lya galaxies, a neutral IGM 
will significantly attenuate their Lya flux. Therefore, the lack of evolution in the Lya galaxy 
luminosity function is used to put an upper limit on the IGM neutral fraction 50% at z ~ 6.5. 

Hainran (2002) gives the radius of an HII region R s in a uniform IGM with neutral fraction 
/hi around a star forming galaxy with star formation rate of M*, 
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A L 


1/3 


40 M 0 yr -1 


10 7 yi 


1/3 


fl/3 f— 1/3 
Je iffl 


1 + Z 

7 


-i 


proper Mpc, 


( 15 ) 


where f* is the duration of the starburst, and f e is the escape rate of Lyman Limit photons. 
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Fig. 12.— Average length of dark gaps with r > 3.5 as a function of IGM neutral fraction, assuming 
that transmitting regions in the quasar spectrum are generated by HII regions surrounding early 
star forming galaxies which follow the Yan & Windhorst (2004b) luminosity function. The solid 
(red) curve is the default case where the duration of star formation in each galaxy is 10 7 yr and the 
escape fraction of ionizing photons is 10%. The two groups of dashed curves above and below the 
default case correspond to the number of escaping photons 3 and 10 times larger or smaller than 
the default case. The dashed line is the observed limit on the average dark gap length at z ~ 6, 
consistent with a neutral fraction <; 0.5. 
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For such an HII region, the GP damping wing optical depth at a distance r from the center 
of the galaxy is approximately (Cen 2003c): 

Ti{T) = L2/hi (16) 

at large distances from the center of the resonant Lya line. In order to preserve a dark gap 
threshold r m i n > 3.5, we require the center of the galaxy HII region to have an optical depth from 
the damping wing, r > 3.5. For a neutral IGM fni = 1, we have R s = 0.34 Mpc, corresponding 
to a star forming rate of M* ~ 70 Mgyr ^ 1 using Eq. (15), assuming that the duration of star 
formation is i* = 10 7 yr and the escape fraction is f e = 0.1 (Steidel et al. 2001). This is a luminous 
galaxy at z ~ 6 : using the calibration of Bunker et al. (2004), it has a zab magnitude of 24.3. 
Such galaxies are extremely rare at z ~ 6 and are not likely to generate transmitting spikes every 
~ 100 comoving Mpc. The fact that the gap length is ~ 100 Mpc therefore already rules out a 
completely neutral IGM. To include the contribution of faint galaxies and place a more stringent 
constraint on the neutral fraction, we assume a luminosity function of star forming galaxies at 
z ~ 6 following Yan & Windhorst (2004b): 

$(Mi3oo) = 0.9214>*10 _0 - 4( “ +1)(M i3oo- M i3oo) exp[-10°' 4(M i3oo- M i3oo)] ) ( 17 ) 

where = 4.55 x 10 - 4 Mpc -3 , a = —1.6 is the faint end slope of the luminosity function, M 1300 
is the absolute AB magnitude at rest-frame wavelength 1300A, and the characteristic luminosity 
Al 1*300 = —21.03. The average length of the dark gap is then the distance between two star forming 
galaxies that have r < 3.5 at the center of their respective HII regions: 

( 1 /A r ) = J 4>(.Mi3oo)7r73 5 (Mi3oo)dMi3oo, (18) 

where r 3 .s(Mi 3 oo) is the radius up to which the damping wing optical depth is still smaller than 
3.5, and the integral is over the luminosity range in which the Lya transition from the galaxy can 
survive the damping wing absorption. We use the calibration of Bunker et al. (2004) to convert 
UV luminosity to star formation rate. 

Figure 12 shows the gap length Ar(r > 3.5) as a function of neutral fraction. In the case of 
f* = 10' yr and f e = 0.1, the average gap length of ~ 80 comoving Mpc at z ~ 6 corresponds 
to a neutral fraction of ~ 0.2. For a population of star forming galaxies with 10 times more 
photon production (oc f e t *), the same gap length corresponds to fni ~ 0 . 6 , while for a population 
with 10 times less photons, the constraint on fni ~ 0.05. The constraint in this simple model, 
however, is only an upper limit on fni for several reasons: ( 1 ) we do not consider the fact that the 
recombination timescale in the early Universe is long compared to the Hubble time, and assume 
that all IGM ionization is produced by HII regions around star forming galaxies. In reality, the 
IGM is probably already largely ionized by previous generations of star forming galaxies; there are 
highly ionized regions without active star forming galaxies in them at the observed epoch (Wyithe 
& Loeb 2005); (2) we do not consider the dumpiness of the IGM. Early galaxies are formed in 
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the highest density, most biased peaks in the density field, and have the highest IGM density 
in the surrounding region. So the constraint here is closer to the neutral fraction in the densest 
part of the universe, which is higher than the average neutral fraction; (3) we do not consider the 
clustering of star forming galaxies at high redshift. If they are strongly clustered, the distance 
between separate clusters of star forming galaxies is longer for the same IGM neutral fraction. A 
lower IGM neutral fraction is needed to explain the observed gap length. Therefore, we regard 
the constraint given in Figure 12 to be an upper limit: the finite distance between dark gaps at 
z ~ 6 suggests that the IGM is not likely to be completely neutral in the regions between HII 
regions ionized by early star forming galaxies. However, a neutral fraction as high as 10-50% is 
still permitted in the optimal case. 


6. HII Regions Around Quasars 


UV photons from a luminous quasar will ionize a HII region (Madau & Rees 2000, Cen & 
Hairnan 2000) in the IGM. The evolution of a cosmological HII region expanding into a neutral 
IGM is described by the equation (e.g., Shapiro & Giroux 1982, Madau, Haardt & Rees 1999): 


n H 




= Nq~ 


nHVi 

tree 


(19) 


where Vj is the HII region volume, and Nq is the number of ionizing photons the central quasar 
emits per unit time. When the quasar lifetime is much less than Hubble time, as expected here, 
the Hubble expansion term is not important. When the HII region is expanding into a partially 
ionized IGM, as in the case discussed in this paper, (1) the left hand side of Eq. (19) becomes 
d(fmn-HVi)/dt, (2) the first term of the right hand side of Eq. (19) becomes Nq + esVi, where 
€b is the emissivity of the diffuse UV background. Assuming photoionization-recombination 
equilibrium in the IGM outside quasar HII region: e# = nHiin e a(T) = (1 — fm) 2 hi/trec (for 
simplicity, we assume a pure hydrogen IGM here), Eq. (19) then changes to: 
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The quasar ionizing photons are only responsible for ionizing a fraction /hi of the IGM, while the 
rest is ionized by the diffuse UV background same as the case outside quasar HII region. The factor 
of (2 — /hi) in Eq. (20) arises from the (1 — /hi) 2 dependence in photoionization-reconrbination 
equilibrium equation. In the case of a mostly ionized IGM, /hi < 1, we have 
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Therefore, for the same quasar lifetime t, the HII region size scales as 
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This scaling applies for both the cases in which quasar lifetime is longer or shorter than the 
recombination time. When recombination is not important (Madau & Rees 2000, Cen & Haiman 
2000, cf. Yu & Lu 2005), Cen & Haiman (2002) show that the HII region has a radius 

R s = 8.0f H )^ 3 (N Q /6.5 x 10 57 s' 1 ) 1/3 (tQ/2 x 10 7 yr) 1/3 [(l + zq)/7] -1 proper Mpc. (23) 

This is the same as Eq. (15), now expressing the source luminosity in terms of the photon 
production rate. The SDSS z ~ 6 quasars would produce HII regions with R s of the order several 
Mpc if the IGM were largely neutral (e.g., Pentericci et al. 2002, White et al. 2003). A number 
of authors (Mesinger & Haiman 2004, Wyithe & Loeb 2004a, Wyithe et al. 2005, Yu & Lu 2005) 
used the observed size of the HII regions around a number of SDSS z ~ 6 quasars to constrain the 
IGM neutral fraction, and suggested that the IGM could have neutral fraction ^ 0.2 outside the 
quasar HII region. Using this size to derive the IGM neutral fraction, however, is complicated by 
a number of factors: 

1. Uncertainties in the quasar lifetime. Current observations suggest that quasar lifetimes lie 
in the range 10 5 - 10 8 years (e.g. Hopkins et al. 2005, Martini & Weinberg 2001). Most 
studies have assumed a lifetime of ~ 10' yrs, similar to the Eddington timescale for black 
hole growth. For a given R s , /hi ~ t.Q when recombination is not important. 

2. Uncertainties in the ionizing photon production rate by the quasar for a given UV luminosity. 

3. Radiative transfer and dumpiness of the IGM within the quasar HII region, which will affect 
the apparent line of sight size of the HII region. 

4. Contribution of ionizing photons from galaxies in the overdense environment of luminous 
quasars. 

5. Large scale clustering of the IGM around luminous quasars at z ~ 6, which are highly 
biased, will affect the estimates of the IGM neutral fraction. 

The last two factors could result in significant overestimation of the neutral fraction. The 
SDSS quasars considered here are very rare objects (spatial density ~ 10 -9 Mpc -3 ) and thus 
are likely to live in highly biased environment in the early Universe. Fan et al. (2001) estimated 
that these quasars correspond to 5-6cr peaks in the density field. Yu & Lu (2005) argued that up 
to ~ 90% of the neutral hydrogen could have already been ionized by the star forming galaxies 
associated with this high density peak. Mesinger & Haiman (2004) showed that by comparing 
high S/N, high resolution spectra with detailed hydrodynamic simulations of quasar HII regions, 
it is possible to reliably estimate the neutral fraction and minimize the degeneracies with quasar 
lifetime and bolometric luminosity along some lines of sight. 

For these reasons, the absolute measurement of the neutral fraction based on HII region size 
could be off by a large factor. Over the narrow redshift range considered here, these systematics 
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are likely to be roughly the same on average for all the quasars. If there is an order of magnitude 
evolution in the IGM ionization, the size of the HII region should show strong evolution, providing 
a reliable relative measurement of neutral fraction. In this section, we study the distribution and 
evolution of HII region sizes around the nineteen z ~ 6 quasars in our sample. 

We first examine the definition and measurement of the HII region size around quasars. The 
physical picture is that quasars are located within an IGM with neutral fraction /hi- As the 
quasar evolves, its ionizing flux produces a surrounding region within which the IGM is completely 
ionized. We have shown above that for a partially ionized IGM, the size of this proximity zone 
grows as ~ Nq 3 t l q 3 . In the case that the recombination timescale is long, this is not a picture of 
a static Stromgren sphere. When the general IGM is already highly ionized (valid at least up to 
z ~ 6, as we have seen), it is not even a standard HII region. It is also very different from the 
classic foreground proximity effect at lower redshift (e.g. Bajtlik, Duncan, & Ostriker 1988, Scott 
et al. 2002). At z < 5, the IGM ionizing background is much higher. At a distance R s from the 
quasar, the UV background itself is much larger than the flux intensity produced by the quasar in 
question. Furthermore, the mean free path of ionizing photons is much larger than R s calculated 
above. In this case, the expanding proximity zone does not have a clear boundary. The proximity 
effect is simply an excess of ionizing background and increased ionization close to the quasar. 

At z ~ 6, with a low UV background, the quasar flux does produce a highly ionized proximity 
zone with a well-defined boundary expanding into the general IGM. However, the boundary of 
this proximity zone can be difficult to measure observationally. It cannot be simply defined as the 
point at which the transmitted flux reaches zero for the following reasons: (a) if the IGM is highly 
ionized, not only does the transmission never reach zero on average, it has also large fluctuations 
(§3). Simply defining the size of the proximity zone as the point at which the flux reaches zero 
will overestimate its size at low redshift. (b) The IGM transmission is characterized by a mixture 
of dark gaps and sharp transmission peaks (§5). Using zero flux as the boundary is affected by the 
resolution and smoothing length used, (c) In the case of a neutral IGM, reaching zero flux does 
not necessarily mean reaching the boundary of the proximity zone: the IGM could still be ionized 
by the quasar to a level much higher than the surrounding material, but with optical depth higher 
than the detection limit (/hi > 1CD 3 ). 

Therefore, instead of using the Stromgren sphere radius R s , we define the size of the 
proximity zone R p as the region where the transmitted flux ratio is above 0.1, when smoothed 
to a resolution of 20A. For the quasars in our sample (z > 5.7), the average IGM transmission 
is smaller than 0.04. Therefore the proximity zone defined this way is not affected by the IGM 
properties outside the quasar HII region, even with the large observed fluctuation in the IGM 
transmission. To first order, at R p , the IGM ionization state is the same for all the quasars, so we 
have: R p oc (Nqtq) 1 ^ (/hi) 1,/3 • The quasars in our sample have very similar intrinsic properties. 
So even though we cannot calibrate this relation on an absolute scale, we can use the evolution of 
R p to measure the relative change in the IGM neutral fraction into which the quasar HII region is 
expanding. 
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Fig. 13.— Size of quasar proximity zone as a function of the redshift of the quasar. The size R p 
is the line of sight distance from the quasar to the point at which the transmitted flux ratio falls 
to 0.1 of the continuum level. The radius also has been scaled to a common absolute magnitude 
(M i 45 o = —27). The linear fit shows a factor of 2.8 decrease in the HII region size around quasars 
from z = 5.7 to 6.4. The error on R p is dominated by the uncertainty in the quasar redshift. 
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The other uncertainty comes from the difficulty in determining the quasar redshift. It is well 
known that there is a systematic velocity offset between high-ionization emission lines such as 
CIV and SilV and the systematic redshift of the quasar host galaxy, as measured by narrow lines 
such as [OIII], or CO molecular lines in the host galaxy. On the other hand, low-ionization lines 
such as Mgll show little offset from the systematic redshift. Richards et al. (2002b) used a large 
sample of SDSS quasars to show that CIV is blueshifted from Mgll by 824 ± 511 km s -1 . At 
z ~ 6 , we have ZMgii = zqiv + (0.019 ± 0.012). Most of the redshifts in Table 1 are determined 
using high-ionization lines such as CIV, SilV or NV. When calculating the size of the proximity 
zone, this shift becomes a major source of systematic error. Five of the quasars in our sample have 
Mgll or CO redshifts: ,11148+5251: zqo = 6.42 (Walter et al. 2003), J1623+3112: ZMgii = 6.22 
(L. Jiang et al. in preparation), J1630+4012: ZMgii = 6.07, J1030+0524: ZMgii = 6.31 (Iwamuro 
et al. 2004), and J0836+0054: ZMgii = 5.82 (Pentericci et al. 2005). We adopt the Mgll or CO 
redshift for these five objects when calculating the size of the proximity zone. For the rest of the 
sample, we assume ZMgll = ^high-ionization + 0.02 in the calculation, and an uncertainty in the 
redshift determination of 0 . 02 . 

Figure 13 shows R p as a function of redshift. We have scaled R p to a common 
absolute magnitude M 1450 = —27 (assumed to be proportional to N q ) as R p (corrected) 

= Rp x 10 °' 4 ( 2 1 +A+ 450 )/ 3 . Correcting for luminosity differences significantly tightens the scatter 
around a linear regression. We also exclude the two BAL quasars as well as SDSS J1335+3533 
(z = 5.94) from this test. The latter object is a lineless quasar presented in Fan et al. (2005), in 
which we used the size of the proximity zone to estimate its redshift. The error bars on R p are 
dominated by the uncertainties in the emission line redshift of the quasars a(z ern ) ~ 0 . 02 , yielding 
an uncertainty of ~ 1.2 Mpc. The solid line in Figure 13 is the best linear fit of R p as a function 
of redshift: 

R p = 6.7 — 8.1 x (z — 6 ) Mpc. (24) 

This linear relation has a Pearson correlation coefficient of 0.72, which is significant at the 99.8% 
level for our sample of 16 quasars. We use this linear fit only to establish the correlation between 
R p and redshift, rather than making any detailed statement about the functional dependence of 
R p on z. From z = 5.7 to 6.4, the average size of the proximity zone R p has decreased from 9.1 
Mpc to 3.5 Mpc, a factor of 2.8. From Figure 13, the scatter of HII region sizes at any given 
redshift is of the order of 2 Mpc, not significantly larger than the size of the error bar. Based 
on Eq. (23), this indicates a small dispersion in intrinsic quasar properties, such as bolometric 
correction, quasar lifetime and relative bias of the quasar environment, at a given redshift. 

Figure 14 shows the mean absorption profiles in the quasar proximity zones as a function of 
redshift. We do not scale the profile according to quasar luminosity in this case. The profile is 
noisy due to the small number of quasars (4-8) in different redshift bins and the presence of strong 
absorption lines within the proximity zones, due to the clumpy IGM close to the quasar. However, 
the evolutionary trend is evident: the z = 5.7 — 5.9 profile is the most extended; the transmission 
declines to ~ 0.1 at ~ 8 — 10 Mpc. The z = 5.9 — 6.1 profile reaches the same level at R p ~ 7 Mpc, 
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Fig. 14.— Average absorption profiles in the proximity zone for quasars at three different redshift 
bins: z = 5.7 — 5.9, z = 5.9 — 6.1 and z > 6.1. The Lya +NV emission profile has been fitted with 
a Gaussian and removed. The absorption profile is significantly narrower at z > 6.1 than at lower 
redshifts. 
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while the profile at z > 6.1 is the narrowest, declining to T = 0.1 after only ~ 3 Mpc. This trend 
is not significantly affected by the choice of binning. Also note that the profiles in the low redshift 
bins continue to decline beyond R p until reaching the IGM level, but the level never reaches zero 
flux. Only in the highest redshift bin does the profile reach zero flux, as the IGM outside the 
quasar HII region has r ^ 5. For z > 6, when the flux outside the HII region does reach zero, the 
proximity zone size R p defined this way is ~ 20% smaller than that using the zero-flux size. 

The significant evolution of the size of quasar HII regions validates the physical picture that 
quasar HII regions are expanding into an increasingly neutral IGM at higher redshift. Since 
Rp oc [(1 + z)f hi] -1//3 (after scaling out the dependence on quasar luminosity), a decrease of HII 
region sizes by a factor of 2.8 from z = 5.7 to 6.4 suggests an increase of the neutral fraction by a 
factor of 14. Note that this relation applies when /hi <C 1. as in the case here; at larger nuetral 
fraction, a detailed solution of Eq. (20) is required. This constraint on the relative increase of the 
neutral fraction does not depend on the model of photoionization and IGM density distribution 
we used in §3. 

Mesinger & Hainran (2004) examined the absorption profile in the line of sight of SDSS 
J1030+0534 (z = 6.28). In this line of sight, the Ly/3 transition reaches zero flux at a redshift 
~ 0.02 lower than that of the Lya transition. They model this feature as the boundary of a 
Stromgren sphere, and suggest that the size of this Stromgren sphere is consistent with a largely 
neutral IGM. From Figure 10, we clearly see the extra transmission in Ly (3 for SDSS J1030+0524. 
However, we do not observe this feature in any other line of sight. 

Our GP optical depth measurements show that at z a f, s ~ 5.7, the volume- and mass-averaged 
neutral fractions are ~ 9.3 x 10" 5 and 2.8 x 10 3 respectively (see Figure 7). Our result here 
implies that the neutral fraction has increased by a factor of 14 at z 6.4, to 1.3 x 10 3 and 
0.04, respectively. Note that the high mass averaged neutral fraction of 0.04 is lower than several 
previous estimates (Mesinger &; Hainran 2004, Wyithe et al. 2005) based on measurements of HII 
region size, but still is in general agreement. Our estimate of neutral fraction evolution is based 
on a simple model in which a uniform UV background ionizes the IGM to a neutral fraction of 
/hi, and the quasar provides the extra UV photons that produce a completely ionized HII region. 
Yu & Lu (2005) describe a more detailed, self-consistent model that includes ionization from 
both stellar sources and quasars. Their model considers in detail the influence of the bias factor 
of quasar environment and quasar lifetime on the size of the quasar HII region. Such detailed 
modeling is beyond the scope of this paper. 


7. Discussion and Summary 

The goal of this paper is to study the evolution of the IGM ionization state near the end of 
reionization using absorption spectra of the highest redshift quasars known. We use three types 
of measurements to study IGM properties. We now discuss their advantages and disadvantages, 
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with an eye towards application to yet higher redshift. 

1. Evolution of Gunn-Peterson optical depth. In §3, we measured the GP optical depth in 
the Lya , Ly f3 and Lyy transitions. This is directly linked to the neutral hydrogen density of the 
IGM. Using a photoionization model, we can compute the underlying UV ionizing background, 
and estimate the neutral fraction and mean free path of UV photons. However, as the optical 
depth increases to r ^ 5 at z > 6, the Lya transition becomes completely saturated. It is 
difficult to place stronger constraints as the measurement begins to be dominated by systematic 
observational errors in the red and near-IR part of the spectrum. The Ly/3 and Lyy transitions 
(§3.2) give more stringent constraints due to their smaller oscillator strengths, but they can 
only improve the limit by a factor of 2 — 4. In addition, one has to correct for the foreground 
absorption from lower order transitions. Direct measurement of the GP optical depth can only 
constrain the neutral fraction of the IGM to lower limits in the range 10 3 to 10 2 . Since at 
high redshift, all IGM transmission is through a small number of transmitting spikes with high 
ionization and/or low density, the effective optical depth is actually not sensitive to the region 
of the IGM where most of the neutral hydrogen lies. Furthermore, the interpretation of optical 
depth depends on the photoionization model and IGM density distribution models used. The 
exact values of derived ionization parameters, such as ionizing background and neutral fraction, 
are therefore model-dependent. But for a narrow redshift range, our calculation shows a robust 
strong evolutionary trend on the ionization parameters. Thus using the evolution of GP optical 
depth is most useful at the low redshift end, and it begins to lose its usefulness when the neutral 
fraction increases to ~ 10~ 3 . 

2. Distribution of Gunn-Peterson troughs and dark gaps. In §4, we discussed using the 
distribution of absorption gaps to characterize the evolution of the IGM. Qualitatively, these 
statistics reveal very similar trends to the GP optical depth measurements. But they also contain 
high order information, related to the clustering of transmission spikes and dark gaps. The 
evolution of dark gaps appears to be the most sensitive quantity tracing the reionization process: 
of the statistics we examine, it shows the most dramatic evolution at z > 6. The distribution 
of dark gap length is readily computable from both observations and cosmological simulations. 
Simulations show that even when the average neutral fraction is high, the dark gap length is 
still finite due to the presence of regions that were ionized by star-forming galaxies. In §5.3, 
we used the length of dark gaps and the observed galaxy luminosity function at z ~ 6 to put 
an independent upper limit on the IGM neutral fraction. Of course, these statistics are not 
straightforward to calculate from analytic or semi-analytic models, and must be derived from 
detailed simulations that include the clustering of the IGM and ionizing sources, and take into 
account spectral resolution and S/N of observed spectra. 

3. HII regions around high-redshift sources. In §6, we measured the HII region sizes around 
high-redshift quasars, and showed that the sizes decrease significantly with redshift, consistent with 
a rapid increase in the IGM neutral fraction. This method provides a powerful and independent 
way of measuring the IGM neutral fraction, and can in principle be applied to regions with high 
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neutral fraction if luminous sources can be discovered. The disadvantage is that the sizes of HII 
regions are quite model dependent: theoretical calculations require not only knowledge of IGM 
properties, but also the intrinsic properties of quasars, including their lifetime and bias relative to 
the density field at high redshift. Further theoretical modelling will improve this constraint. As 
these quantities are unlikely to change over the redshift range we probe, we used the HII region 
sizes as a robust measure of the fractional change in the neutral fraction over this redshift range. 

In §1, we posed three questions regarding the evolution of the IGM and the end of reionization. 
We now present our major conclusions based on our analysis of this sample of nineteen z > 5.7 
quasars. 

1. There is a strong evolution of the IGM ionization state at high redshift, in terms of the 
effective GP optical depth, the derived UV background, the neutral fraction, the mean free 
path of UV photons, and the extent of dark gaps. This evolution accelerates at z > 5.7. The 
GP optical depth evolution changes from Tq® ~ (1 + z ) 4 ' 3 to (1 + z)^ 11 , indicating a lower 
limit in the increase in the volume-averaged neutral fraction of a factor of 7 between z = 5.5 
and 6.2. 

2. Accompanying the increased average neutrality of the IGM is a rapid increase in the relative 
dispersion of IGM properties along different lines of sight. We find that this increased 
dispersion implies fluctuations by a factor of ;> 4 in the UV background at z > 6. This is not 
surprising, as the mean free path of UV photons at z ~ 6 is comparable to the correlation 
length of star forming galaxies which likely provide most of the photons that ionized the 
Universe. 

3. The most effective measure of the IGM evolution at high redshift is the distribution of dark 
absorption gaps in the IGM. The average length of dark gaps shows a dramatic increase 
at z ~ 6. This quantity shows the largest line of sight variations: dark gaps as long as 50 
comoving Mpc appear by z ~ 5.5, while some lines of sight at z > 6 are still somewhat 
transparent. 

4. Using the evolution of HII region sizes around the quasars, we find the neutral fraction of 
the IGM has increased by ~ 14 between z = 5.5 and 6.4, consistent with our GP optical 
depth estimates. This implies a mass-averaged neutral fraction as high as ~ 4% in regions 
around luminous quasars at z ~ 6.4. 

5. There is no compelling evidence that the IGM at z ~ 6.4 is largely neutral. On the contrary, 
the finite length of dark gaps in the spectra implies an upper limit on the neutral fraction 
less than 50%, likely even lower. 

Figure 15 summarizes our constraints on the nrass-averaged neutral fraction using quasar 
spectra presented in this paper. The analysis here confirms the main result of Paper I based on 
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Fig. 15.— The evolution of the mass-averaged neutral fraction using constraints from GP optical 
depth measurements, HII region size evolution, and distribution of dark absorption gaps. 
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only four quasars: the detections of deep GP troughs indicate an accelerated rate of IGM evolution 
at z > 5.7, consistent with the IGM transition at the end of the overlapping stage of reionization. 
However, the results do not indicate that the IGM has achieved a high level of neutrality at z ~ 6. 
The large dispersion of the IGM properties along different lines of sight strongly suggests that 
the reionization process is complex and not likely a uniform phase transition over a very narrow 
redshift range. 

If luminous sources at z ~ 6 — 10 can be discovered by JWST and next generation wide-held 
IR surveys, statistics of dark gaps in their absorption spectra and measurements of HII region sizes 
surrounding these sources will provide the most useful probe of the history of reionization using 
high-redshift discrete sources, and complement constraints from CMB polarization measurements 
(e.g. PLANCK) and 21cm experiments (e.g. LOFAR, MWA, PAST and SKA). 
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Table 1. Spectroscopic Observations of the z ~ 6 Quasar Sample 


Quasar (SDSS) 

redshift 

A7l450 

Telescope 1 

Exposure 

(hours) 

Ref 2 

J114816.64+525150.3 

6.42 

-27.8 

Keck 

11.5 

4,5 

J103027.10+052455.0 

6.28 

-27.2 

Keck 

10.3 

2,5 

J162331.81+311200.5 

6.22 

-26.6 

Keck 

1.0 

6,8 

J104845.05+463718.3 3 

6.20 

-27.6 

Keck 

4.5 

4,8 

J125051.93+313021.9 

6.13 

-27.1 

Keck 

1.0 

7,7 

J160253.98+422824.9 4 

6.07 

-26.8 

MMT 

1.3 

6,6 

J 163033.90+401209.6 

6.05 

-26.1 

HET 

1.0 

4,4 

J113717.73+354956.9 

6.01 

-27.1 

Keck 

0.7 

7,7 

J081827.40+172251.8 

6.00 

-27.4 

Keck 

0.3 

7,7 

J 130608.26+035626.3 

5.99 

-26.9 

Keck 

0.2 

2,3 

J133550.81+353315.8 

5.95 

-26.8 

MMT 

2.0 

7,7 

J141111.29+121737.4 

5.93 

-26.8 

Keck 

1.0 

6,8 

J084035.09+562419.9 

5.85 

-26.6 

MMT 

1.0 

7,7 

J000552.34—000655.8 

5.85 

-26.5 

Keck 

0.3 

6,8 

J143611.74+500706.9 

5.83 

-26.3 

MMT 

1.0 

7,7 

J083643.85+005453.3 

5.82 

-27.8 

Keck 

0.3 

2,3 

J000239.39+255034.8 

5.80 

-27.7 

MMT 

1.0 

6,6 

J092721.82+200123.7 

5.79 

-26.8 

KP 

3.7 

7,7 

J104433.04-012502.2 3 

5.74 

-27.5 

Keck 

0.7 

1,1 


1. Telescopes and instrument: Keck: Keck/ESI; MMT: MMT/Red Channel; KP: Kitt Peak 
4-meter/MARS 

2. References of discovery spectrum and spectrum presented here; (1). Fan et al. 2000, (2). 
Fan et al. 2001, (3) Becker et al. 2001, (4). Fan et al. 2003, (5). White et al. 2003, (6). Fan et al. 
2004, (7). Fan et al. 2005, (8). This paper 

3. BAL quasar 

4. This object was listed as J160254.18+422822.9 in Fan et al. (2004); the position given here 
reflects improved astrometry. 
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Table 2. Lya Transmitted Flux Ratio 


Quasar 

4m 

Zabs 

T 

J0002+2550 

5.80 

5.58 

0.0170 ± 0.0062 



5.43 

0.0573 ± 0.0066 



5.28 

0.0205 ± 0.0045 



5.13 

0.1243 ± 0.0050 



4.98 

0.1002 ± 0.0050 

J0005—0006 

5.85 

5.64 

0.0823 ± 0.0069 



5.49 

0.0718 ± 0.0070 



5.34 

0.0961 ± 0.0066 



5.19 

0.0578 ± 0.0063 



5.04 

0.1567 ± 0.0073 

J0818+1722 

6.00 

5.81 

0.0216 ± 0.0053 



5.66 

0.0440 ± 0.0040 



5.51 

0.0984 ± 0.0043 



5.36 

0.1192 ± 0.0039 



5.21 

0.0884 ± 0.0039 



5.06 

0.1285 ± 0.0042 

J0836+0054 

5.82 

5.52 

0.0907 ± 0.0011 



5.37 

0.0348 ± 0.0009 



5.22 

0.0606 ± 0.0009 



5.07 

0.0751 ± 0.0011 



4.92 

0.1276 ± 0.0011 

J0840+5624 

5.85 

5.66 

0.0883 ± 0.0176 



5.51 

0.1127 ± 0.0260 



5.36 

0.1661 ± 0.0202 



5.21 

0.1191 ± 0.0167 



5.06 

0.1765 ± 0.0190 

J0927+2001 

5.79 

5.61 

0.0884 ± 0.0136 



5.46 

0.1041 ± 0.0163 



5.31 

0.0596 ± 0.0104 



5.16 

0.1165 ± 0.0105 



5.01 

0.1268 ± 0.0117 

J1030+0524 

6.28 

6.10 

0.0012 ± 0.0010 



5.95 

0.0060 ± 0.0010 



5.80 

0.0260 ± 0.0012 



5.65 

0.0462 ± 0.0009 



5.50 

0.0661 ± 0.0009 



5.35 

0.1147 ± 0.0008 
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Table 2. Lya Transmitted Flux Ratio (continued) 


Quasar 

4m 

Zabs 

T 

J1044—0125 

5.74 

5.55 

0.0686 ± 0.0022 



5.40 

0.0520 ± 0.0020 



5.25 

0.0427 ± 0.0019 



5.10 

0.0898 ± 0.0020 



4.95 

0.1139 ± 0.0022 

J1048+4637 

6.20 

5.68 

0.0117 ± 0.0011 



5.53 

0.0519 ± 0.0012 



5.38 

0.0736 ± 0.0011 

J1137+3549 

6.01 

5.83 

0.0116 ± 0.0029 



5.68 

0.1010 ± 0.0024 



5.53 

0.0742 ± 0.0026 



5.38 

0.1341 ± 0.0022 



5.23 

0.1323 ± 0.0023 



5.08 

0.0530 ± 0.0025 

J1148+5251 

6.42 

6.25 

0.0015 ± 0.0005 



6.10 

0.0051 ± 0.0005 



5.95 

0.0038 ± 0.0005 



5.80 

0.0186 ± 0.0006 



5.65 

0.0433 ± 0.0005 



5.50 

0.0278 ± 0.0005 

J1250+3130 

6.13 

5.90 

0.0108 ± 0.0033 



5.75 

0.0055 ± 0.0030 



5.60 

0.0248 ± 0.0026 



5.45 

0.0077 ± 0.0026 



5.30 

0.0776 ± 0.0024 

J1306+0356 

5.99 

5.77 

0.0645 ± 0.0020 



5.62 

0.0690 ± 0.0016 



5.47 

0.0991 ± 0.0015 



5.32 

0.0864 ± 0.0014 



5.17 

0.1156 ± 0.0013 
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Table 2. Ly a Transmitted Flux Ratio (continued) 


Quasar 

Zem 

Zabs 

T 

J1335+3533 

5.94 

5.73 

0.0224 ± 0.0059 



5.58 

0.0445 ± 0.0074 



5.43 

0.1215 ± 0.0083 



5.28 

0.1217 ± 0.0058 



5.13 

0.1293 ± 0.0064 

J1411+1217 

5.93 

5.71 

0.0322 ± 0.0033 



5.56 

0.0665 ± 0.0031 



5.41 

0.0858 ± 0.0029 



5.26 

0.0690 ± 0.0028 



5.11 

0.1650 ± 0.0030 

J1436+5007 

5.83 

5.66 

0.0714 ± 0.0217 



5.51 

0.0775 ± 0.0315 



5.36 

0.0895 ± 0.0239 



5.21 

0.1292 ± 0.0199 



5.06 

0.1509 ± 0.0224 

J1602+4228 

6.07 

5.85 

0.0687 ± 0.0057 



5.70 

0.0729 ± 0.0044 



5.55 

0.0795 ± 0.0058 



5.40 

0.0802 ± 0.0055 



5.25 

0.0934 ± 0.0036 

J1623+3112 

6.22 

6.08 

-0.0071 ± 0.0020 



5.93 

0.0125 ± 0.0022 



5.78 

0.0071 ± 0.0024 



5.63 

0.0402 ± 0.0018 



5.48 

0.0407 ± 0.0017 



5.33 

0.0546 ± 0.0016 

J1630+4012 

6.05 

5.77 

-0.0165 ± 0.0342 



5.62 

0.0495 ± 0.0323 



5.47 

0.1015 ± 0.0353 



5.32 

0.1376 ± 0.0296 



5.17 

0.0869 ± 0.0219 











Table 3. Ly/3 Transmitted Flux Ratio 


Quasar 

4m 

Zabs 


T 


J0002+2550 

5.80 

5.58 

0.0557 

± 

0.0039 

J0005—0006 

5.85 

5.64 

0.0786 

± 

0.0066 

J0818+1722 

6.00 

5.81 

0.0341 

± 

0.0035 

J0840+5624 

5.85 

5.66 

0.0898 

± 

0.0127 

J0927+2001 

5.79 

5.61 

0.0715 

± 

0.0088 

J1030+0524 

6.28 

6.10 

-0.0005 

± 

0.0010 

J1044—0125 

5.74 

5.55 

0.0584 

± 

0.0019 

J1137+3549 

6.01 

5.83 

0.0461 

± 

0.0021 

J1148+5251 

6.42 

6.25 

0.0038 

± 

0.0005 



6.10 

0.0165 

± 

0.0005 

J1250+3130 

6.13 

5.90 

0.0128 

± 

0.0024 

J1306+0356 

5.99 

5.77 

0.0458 

± 

0.0014 

J1335+3533 

5.94 

5.73 

0.0406 

± 

0.0058 

J1411+1217 

5.93 

5.71 

0.0300 

± 

0.0028 

J1436+5007 

5.83 

5.66 

0.0584 

± 

0.0180 

J1602+4228 

6.07 

5.85 

0.0607 

± 

0.0035 

J1623+3112 

6.22 

6.08 

-0.0036 

± 

0.0019 



5.93 

0.0305 

± 

0.0017 

J1630+4012 

6.05 

5.77 

0.0499 

± 

0.0194 


This is the total observed transmission, 
without correcting for the foreground Lya 
absorption. 
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Table 4. Ly 7 GP Optical Depth Measurements 


object 

Zem 

redshift 

range 

r G p using Ly 7 

J1030+0524 

6.28 

6.11 - 

6.17 

> 14.2 1 

J1148+5251 

6.42 

6.25 - 

6.32 

14.7 ± 1.1 

J1623+3112 

6.22 

6.05 - 

6.15 

> 12.2 


1. 2-cr lower limit 


Table 5. Evolution of Average GP Optical Depth and Dispersion 


redshift range 

(tgp) 

<7 (r G p) 

4.90 - 5.15 

2.1 

0.3 

5.15 - 5.35 

2.5 

0.5 

5.35 - 5.55 

2.6 

0.6 

5.55 - 5.75 

3.2 

0.8 

5.75 - 5.95 1 

4.0 

0.8 

5.95 - 6.25 1 

7.1 

2.1 


1. The redshift bin includes 
complete GP troughs, so the value 
is a lower limit 







Table 6. 


Dark Troughs with A z >0.1 in the Highest-redshift Quasars 


object 

4m 

redshift range 

TCP 

using Lya 

PGP 

using Ly/3 

J1148+5251 

6.42 

6.10 - 

- 6.32 

6.2T0.2 

7.7T0.2 

J1148+5251 

6.42 

5.91 - 

- 6.07 

6.1T0.2 

5.5T0.1 

J1030+0524 

6.28 

5.98 - 

- 6.17 

> 6.3 1 

> 10.2 

J1623+3112 

6.22 

5.96 - 

- 6.13 

> 5.4 

> 8.2 

J1048+4637 

6.20 

5.63 - 

- 5.75 

> 5.4 

> 5.7 

J1250+3130 

6.13 

5.69 - 

- 5.85 

> 4.6 

6.5T0.6 


1. 2-cr lower limit 







